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ABSTRACT 


A  new  radio  direction  finding  system  has  been  proposed 
which  offers  advantages  over  current  systems  in  the  areas 
of  bearing  accuracy  and  simplicity  of  antenna  design 
requirements.   This  system  phase  codes  the  outputs  of  a 
linear  array  of  omnidirectional  antennas,  and  then  combines 
these  time-separated  signals  so  that  the  sequences  produced 
uniquely  determine  the  signal  arrival  angle.   A  correlation 
process  using  matched  filters  is  performed  on  the  resulting 
sequences  in  order  to  determine  the  corresponding  arrival 
angle. 

Resolution  of  the  inherent  ambiguities  of  the  system, 
due  to  elevation  angle  and  adjacent  quadrant  azimuth  compo- 
nents is  demonstrated  using  a  two-array  technique.   A 
computer  simulation  of  the  phase  coding/decoding  process, 
using  an  idealized  signal  model  is  used  to  evaluate  the 
performance  of  the  system  for  various  input  noise  conditions, 
array  lengths,  and  input  signal  parameters. 
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I.   INTRODUCTION 

A.   STATEMENT  OP  THE  PROBLEM 

In  order  to  determine  the  geographical  location  of  a 
radio  frequency-emitter,  its  emissions  must  be  detected  and 
processed  simultaneously  by  two  or  more  radio  direction 
finding  (RDF)  systems.   The  unambiguous  determination  of 
the  angle  of  arrival  (AOA)  of  RF  energy  received  from  a 
target  transmitter  enables  each  RDF  system  to  generate  a 
line-of-bearing  along  the  Poynting  vector  associated  with 
the  transmitted  energy.   Ideally,  the  lines-of-bearing  can 
then  be  plotted  on  a  chart  containing  the  known  locations 
of  the  RDF  sites,  and  the  point  of  intersection,  or  the 
area  enclosed  by  three  or  more  lines-of-bearing,  will 
indicate  the  approximate  transmitter  location. 

The  accuracy  with  which  a  target  transmitter  can  be 
located  by  such  an  operation  depends  upon  many  factors. 
The  number  and  precision  of  the  individual  AOA  measurements, 
and  the  relative  positions  of  the  transmitter  and  RDF  sites 
are  probably  the  most  important;  and  each  of  these  factors 
is,  in  turn,  dependent  upon  numerous  other  considerations, 
both  technical  and  non-technical.   The  point  of  intersection 
of  two  near-parallel  lines-of-bearing  (LOB)  will  vary 
significantly  with  relatively  small  standard  deviation 
errors  in  either  LOB,  while  the  intersection  of  two  near- 
perpendicular  lines-of-bearing  will  be  influenced  only 
slightly  by  the  same  errors.   It  is  obvious  that  additional 


AOA  data,  providing  a  wide  dispersion  of  available  lines-of- 
bearing,  can  greatly  enhance  emitter  location  accuracy  over 
a  large  geographical  area. 

Increasing  the  number  and  geographical  dispersion  of 
radio  direction  finding  sites  has  proven  very  effective  in 
the  improvement  of  emitter  location  in  the  HF  communications 
spectrum  (2  -  30  MHz).   The  problem  of  accurate  AOA  measure- 
ments is  particularly  acute  in  this  frequency  range  since 
most  signals  of  interest  are  received  via  ionospheric 
propagation.   The  long  propagation  paths  involved,  and  the 
non-deterministic  nature  of  ionospheric  propagation,  have 
necessitated  the  use  of  probabilistic  techniques  applied  to 
multiple  RDF  measurements.   The  availability  of  multiple 
LOB  inputs  from  RDF  sites  in  a  coordinated  net  allows  a 
central  processor  to  sort  through  the  inputs  using  an 
optimization  scheme,  and  to  choose  those  inputs  which  yield 
a  most  probable  emitter  location  with  a  minimum  area  of 
ambiguity. 

From  a  technical  standpoint,  the  greater  the  number  and 
dispersion  of  RDF  sites,  the  larger  the  statistical  data 
base  from  which  decisions  can  be  made.   However,  from 
economic  and  political  standpoints,  the  number  and  placement 
of  complex,  large-array  RDF  systems  for  operation  in  the  HF 
range  are  very  limited.   The  development  of  a  simpler, 
more  compact  HFDF  system  would  allow  the  strategic  placement 
of  additional  sites  in  those  locations  where  more  complex 
facilities  do  not  exist,  or  are  impractical. 


B.  STATEMENT  OF  PURPOSE 

The  purpose  of  this  paper  Is  to  investigate  the  feasibil- 
ity of  a  new  and  unique  radio  direction  finding  concept,  and 
to  determine  the  applicability  of  this  concept  to  a  realizable 
system  for  operation  in  the  HF  communications  range.   The 
fundamental  concept  was  originally  presented  in  the  Masters 
thesis  of  R.  C.  Todaro  [1],   The  results  of  this  original 
work  indicated  that  the  new  system  concept  offered  definite 
advantages  over  current  RDF  systems  in  the  areas  of  bearing 
accuracy  and  simplicity  of  antenna  design. 

The  specific  areas  of  investigation  described  in  this 
paper  are:   (1)  a  discussion  and  analysis  of  the  ambiguities 
inherent  to  the  system,  (2)  a  qualitative  analysis  of  system 
performance  in  the  presence  of  background  noise;  and, 
(3)  physical  scaling  of  the  system  for  operation  in  the  HF 
frequency  range. 

The  remainder  of  this  chapter  contains  a  brief  descrip- 
tion of  a  generalized  RDF  system,  and  a  comparison  of  the 
RDF  systems  currently  in  use.   Chapter  Two  contains  a 
description  of  the  new  system,  and  a  discussion  of  geometric 
ambiguities.   Chapter  Three  describes  the  computer  model  of 
the  system  and  the  results  obtained  from  the  incorporation 
of  background  noise  and  scaling  parameters. 

C.  GENERALIZED  RDF  SYSTEM 

Before  examining  any  radio  direction  finding  system  for 
its  individual  characteristics,  it  is  useful  to  consider  the 


elements  basic  to  any  system  and  understand  how  they  can 
affect  its  performance.   A  block  diagram  of  a  generalized 
RDF  system  is  shown  in  Figure  1-1. 
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Figure  1-1.   A  GENERALIZED  RDF  SYSTEM 

1.   Signal  Source 

The  signal  source  has  only  three  properties  which 
will  affect  the  operation  of  an  RDF  system.   These  are  its 
placement  on  the  earth,  its  transmission  frequency,  and  its 
type  of  modulation.   Its  location  on  the  surface  of  the 
earth  will  determine  the  great  circle  bearing  between  the 
transmitter  and  the  RDF  location,  but  except  in  very  simple 
cases  this  is  insufficient  information  with  which  to  predict 
the  line  of  bearing  readout  from  the  RDF  system.   The  trans- 
mission frequency  will  affect  the  results  obtained  by 
determining,  in  conjunction  with  the  transmission  medium, 
the  mode  by  which  the  signal  is  sent  from  the  source  to  the 
RDF  site.   The  type  of  modulation  will  have  an  effect  on 
those  systems  which  compare  arrival  times  of  specific 
modulation  events,  or  those  which  integrate  multiple  antenna 
sweeps  over  each  line  of  bearing. 
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2 .   Propagation  Medium 

It  is  possible  for  the  propagation  medium  to  have 
little  effect  on  the  angle  of  arrival  of  the  received  signal, 
but  this  is  not  normally  the  case.   If  the  signal  arrives  at 
the  RDF  site  by  ground  wave  propagation,  then  results  are 
relatively  simple  and  easily  predicted.   Most  of  the  signals 
of  interest  arrive  via  the  ionoxphere,  however,  so  its 
effects  must  be  taken  into  account.   There  has  been  a  great 
deal  of  work  done  to  investigate  the  character  of  the 
ionosphere,  and  much  more  remains  before  accurate  predictions 
of  the  transmission  path  of  a  received  signal  can  be  obtained 
A  large  part  of  the  difficulty  lies  in  the  complexity  of  the 
equipment  necessary  for  measurements  on  the  ionosphere,  and 
the  large  amount  of  time  necessary  to  complete  a  set  of 
measurements.   Since  the  characteristics  of  the  ionosphere 
depend  on  the  sunspot  cycle  as  well  as  seasonal  and  diurnal 
changes,  a  complete  set  of  measurements  on  a  given  path 
could  not  be  obtained  in  less  than  eleven  years. 

If  accurate  transmission  path  prediction  were 
available  it  is  conceivable  that  the  emitter  location 
process  could  be  performed  by  a  single-site  RDF  system. 
That  is,  azimuth  and  elevation  AOA  measurements  could  be 
used  in  conjunction  with  path  prediction  to  compute  the 
point  of  origin  of  received  signals.   Such  a  capability 
would  greatly  enhance  the  emitter  location  capability  of  a 
coordinated  HFDF  net. 
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3.   Signal  Receptor  and  Encoder 

This  is  the  first  section  of  an  RDF  system  over 
which  some  design  control  can  be  exercised,  and  the  choice 
of  this  block  will  determine  to  a  great  extent  the  accuracy 
and  speed  with  which  a  line-of-bearing  can  be  generated  from 
the  received  signal.   The  basic  purpose  of  this  unit  is  to 
sense  the  incoming  signal,  combine  the  proper  components  to 
give  an  output  which  can  be  uniquely  interpreted  as  an  angle 
of  arrival,  and  present  these  components  to  a  code  interpreter 
for  transformation  to  a  form  suitable  for  display.   Some  of 
the  common  types  of  receptors  and  encoders  will  be  discussed 
in  a  later  section  describing  the  RDF  systems  currently  in 
use. 

H .   Code  Interpreter  and  Bearing  Output  Unit 

In  most  RDF  systems  the  functions  of  this  section 
are  performed  jointly  by  an  electronic  indicator  sub-system 
and  the  human  operator.   Ideally,  the  code  interpreter 
should  be  carefully  matched  to  the  characteristics  of  the 
coder  so  the  desired  output  information  can  be  obtained  in 
the  shortest  time  and  with  the  highest  possible  accuracy. 
There  is,  of  course,  no  way  in  which  the  equipment  associated 
with  these  functions  can  increase  the  fundamental  accuracy 
of  the  system  to  a  greater  degree  than  is  allowed  by  the 
components  of  the  receptor  and  encoder,  but  it  should  not 
degrade  their  performance.   An  exact  match  occurs  between 
the  signal  encoder  section  and  the  code  interpreter  of  the 
proposed  system  through  the  use  of  a  matched  filter  technique. 
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That  is,  the  encoding  section  generates  a  phase  coded  sequence 
which  is  uniquely  determined  by  the  angle  of  arrival  of  the 
incoming  signal.   The  inverse  filter  in  the  code  interpreter 
which  provides  an  exact  match  for  this  sequence  will  produce 
a  unique  autocorrelation  (or  cross  correlation)  output. 
This  process  is  described  in  more  detail  in  the  following 
chapter. 

D.   CURRENT  RDF  SYSTEMS 

There  are  four  basic  system  types  which  are  commonly 
used  to  determine  the  AOA  of  transmitted  signals.   The 
choice  of  which  system  to  use  in  any  particular  application 
is  determined  primarily  by  the  frequency  range  of  interest, 
the  degree  of  accuracy  desired,  and  the  type  of  signal 
modulation  expected.   These  four  basic  types  are  as  follows: 

1)  Relative  amplitude 

2)  Relative  phase 

3)  Time  of  Arrival  (TOA) 

4)  Doppler 

1.   Relative  Amplitude  Methods 

By  far  the  most  commonly  used  RDF  systems,  these 
employ  directional  antennas  whose  patterns  may  be  mechan- 
ically or  electronically  rotated.  The  line  of  bearing  is 
then  determined  by  observing  the  angular  position  of  the 
antenna  (or  its  beam  pattern)  when  the  system  response  is 
either  at  a  maximum  or  at  a  minimum  value.  Those  systems 
which  employ  physically  small  antennas  (e.g.,  magnetic-field 
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sensitive  loops)  mechanically  rotate  the  antenna  while 
observing  the  system  output. 

A  loop  antenna  which  is  very  much  smaller  than  a 
wavelength  has  two  sharp  nulls  in  its  pattern,  and  these 
nulls  lie  in  a  plane  normal  to  the  plane  of  the  loop. 
Rotating  an  antenna  of  this  type  and  noting  the  positions 
of  the  nulls  will  give  an  indication  of  the  bearing  of  the 
incoming  signal.   The  bidirectional  ambiguity  can  be  resolved 
by  addition  of  the  signal  picked  up  by  a  separate  omni- 
directional antenna  located  closely  to  the  loop,  but  not 
so  close  as  to  affect  the  pattern  of  the  loop.   Addition 
of  the  antenna  outputs  will  produce  a  cardioid  pattern 
indicating  the  proper  null  to  use  for  the  bearing. 

The  signal  processing  portion  of  this  system  (code 
interpreter  and  bearing  output)  consists  very  simply  of  a 
standard  communications  receiver  whose  output  feeds  the 
deflection  circuits  of  a  circularly  swept  CRT  display. 
Antenna  position  is  synchronized  to  the  CRT  sweep  to 
establish  a  common  rotation  rate  and  a  fixed  reference 
bearing. 

The  primary  disadvantage  of  this  system  is  the 
relatively  low  overall  gain  of  the  antenna,  and  the  fact 
that  the  distinct  feature  of  the  pattern  (null)  occurs 
where  the  gain  is  at  a  minimum.   For  these  reasons  the  loop 
antenna  is  ineffective  in  a  low  signal  level/high  noise 
environment.   In  addition,  the  loop  antenna  is  subject  to 
polarization  errors.   For  ionospherically  propagated  waves, 


which  have  nearly  random  polarization,  these  errors  may  be 
as  large  as  ±35°  making  the  loop  antenna  worthless  for  long 
range  DP  in  the  high  frequency  band  [2]. 

Instead  of  using  a  rotating  element,  two  or  more 
fixed  elements  may  be  installed  with  their  respective  nulls 
symmetrically  spaced.   These  individual  elements  may  be 
loops,  pairs  of  monopoles  (thereby  minimizing  the  polariza- 
tion error  mentioned  above),  or  any  structure  which  will 
produce  a  "figure-eight"  pattern  with  distinct  nulls.   The 
outputs  from  these  elements  are  combined  in  a  special  motor/- 
transformer  in  which  they  are  connected  to  symmetrically 
spaced  stator  windings.   A  rotating  pickup  coil  (rotor)  acts 
as  the  secondary  winding,  and  its  output  is  fed  to  the  RDF 
receiver.   The  amplitude  output  of  the  receiver  can  then 
be  compared  with  the  mechanical  position  of  the  rotor  of 
the  transformer,  or  goniometer,  and  the  bearing  of  the  signal 
obtained.   Bidirectional  ambiguities  can  be  resolved  by 
forming  a  cardioid  pattern  as  in  the  rotating  loop  system. 

While  the  figure  eight  patterns  of  the  individual 
elements  have  relatively  broad  beamwidths,  the  combining 
action  of  the  goniometer  produces  a  sharper  "propellor" 
pattern  with  a  beamwidth  inversely  proportional  to  the 
number  of  individual  element  pairs.   This  allows  the  use  of 
the  maximum  gain  portion  of  the  array  pattern  for  determin- 
ation of  AOA.   However,  in  order  to  achieve  good  figure 
eight  patterns,  small  aperture  elements  (much  less  than  a 
wavelength)  must  be  employed.   Therefore  the  array  still 
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has  a  low  overall  gain  characteristic.   In  addition,  the 
small  aperture  requirement  severely  limits  the  number  of 
non-interacting  antenna  elements  which  may  be  used  to  com- 
prise the  array.   This  limits  the  factor  by  which  the  overall 
antenna  pattern  beam  width  may  be  reduced,  and,  in  turn, 
limits  the  angular  resolving  capability  of  the  array. 

In  order  to  improve  the  gain  and  bearing  resolving 
capabilities  of  relative  amplitude  RDF  systems,  large 
aperture  arrays  have  been  developed.   These  arrays  are 
usually  composed  of  a  large  number  of  individual  antennas 
located  on  the  periphery  of  a  symmetrical  structure,  ordi- 
narily a  circle,  and  connected  to  a  goniometer  which  feeds 
a  single  channel  receiver.   The  code  interpreter/bearing 
output  section  of  this  system  is  similar  in  design  and 
operation  to  that  associated  with  the  simple  spinning  loop 
system. 

The  large  aperture  RDF  arrays  that  are  most  commonly 
in  use  today  are  of  the  Wullenweber  design.   This  is  a 
multiband  circular  array,  each  section  of  which  contains 
an  independent  concentric  array  of  vertical  monopoles,  and 
a  circular  reflecting  screen.   Two  or  three  concentric 
sections,  each  with  different  numbers  and  sizes  of  active 
monopole  elements,  are  designed  to  cover  overlapping  portions 
of  the  HF  communications  spectrum.   The  outputs  from  each 
of  the  individual  antenna  elements  are  fed  to  a  capacitive 
coupled  goniometer,  which  utilizes  discrete  rotor  and  stator 
pickup  plates  instead  of  inductively  linked  coil  windings. 
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The  use  of  capacitive  coupling,  required  by  the  large 
number  of  individual  antenna  elements,  enables  the  code 
interpretor/bearing  output  section  to  take  full  advantage 
of  the  fine-grain  resolving  capability  of  the  antenna  array. 
The  goniometer  forms  a  very  narrow  array  beam  from  a  set  of 
40  pickup  plates  located  on  the  rotor,  with  these  40  plates 
feeding  delay  lines  that  act  to  compensate  for  the  physical 
structure  of  the  array  and  transform  it  into  an  equivalent 
linear  (broadside)  array. 

The  Wullenweber  can  thus  produce  a  very  narrow 
rotatable  beam  providing  resolving  power  of  approximately 
3°.   However,  in  order  to  achieve  this  resolving  power,  up 
to  120  separate  antenna  elements  are  required  for  each 
frequency  range,  and  the  physical  dimensions  of  the  array 
commonly  exceed  1000  feet  in  diameter  and  200  feet  in 
height.   Therefore,  the  array  and  its  associated  transmission 
line  system  is  a  very  complex  and  expensive  structure. 

All  relative  amplitude  systems,  regardless  of  array 
size,  have  the  disadvantage  that  they  are  not  instantaneous. 
That  is,  the  amount  of  time  the  narrow  array  beam  spends  at 
each  azimuth  angle  is  a  small  fraction  of  its  total  rota- 
tional period.   This  can  produce  erroneous  and  misleading 
results  with  signals  of  short  duration  such  as  on/off  keyed 
pulse  transmissions.   The  scan  rate  may  be  increased  to 
overcome  these  errors,  but  a  corresponding  increase  in  the 
RDP  receiver  bandwidth  is  required  in  order  to  recover  the 
rapidly  sampled  signals.   Therefore,  the  maximum  scan  rate 
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is  limited  by  the  receiver  bandwidth  necessary  for  good 
co-channel  discrimination  (selectivity). 
2.   Relative  Phase  Methods 

These  systems  utilize  the  measured  signal  phase 
difference  between  omnidirectional  antennas  to  determine 
angle  of  arrival.   This  phase  difference  is  given  by  the 
equation: 


<j>  =  2tt-7-cos  0    (radians) 


where, 

d  =  antenna  separation  along  the  baseline 

X  =  wavelength  of  incident  wavefront 

6  =  angle  of  arrival,  measured  from  the  baseline. 

It  is  obvious  that,  assuming  a  given  accuracy  in 
phase  measurements,  the  accuracy  with  which  the  AOA  can  be 
determined  (i.e.,  the  resolving  power  of  the  system)  is 
greater  for  larger  values  of  the  ratio  ■*-.      When  r-  is  greater 
than  unity,  however,  there  is  a  possible  ambiguity  of  2tt, 
or  a  multiple  of  2it,  in  the  measured  phase  angle,  and  there 
will  not  be  a  unique  solution  for  the  AOA. 

Providing  that  -r-  is  not  too  large  however,  the 
alternative  solutions  are  well  separated,  and  it  is  gener- 
ally possible  from  other  considerations  to  select  the  most 
probable  [3].   For  example,  the  approximate  wavelength,  or 
range  of  wavelengths,  of  the  incident  signals  can  be 
determined  from  the  tuned  frequency  of  the  RDF  receivers, 
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le  the  separation  (d)  is  a  fixed  constant.   Simple  logic 
cuitry  can  be  employed  to  associate  a  range  of  values  of 

phase  shift  with  corresponding  frequency  bands  of  the 

receivers . 

An  additional  ambiguity  results  when  the  incident 
,nal  wavefront  has  a  non-zero  elevation  angle  component 

to  ionospheric  propagation.   The  phase  difference 
ation  will  contain  an  additional  cosine  factor  of  the 
vation  AOA,  and  the  one-equation/one-unknown  situation 
1  no  longer  apply.   The  requirement  for  an  additional 
sured  parameter  to  uniquely  determine  both  elevation  and 
muth  angles  of  arrival  necessitates  the  use  of  a  second 
n-collinear)  antenna  pair.   The  second  antenna  pair  is 
ally  placed  orthogonally  to  the  first,  resulting  in  a 
xially  symmetric  array.   Since  the  elevation  angle 
ponent  must  be  considered  in  the  proposed  new  system, 
s  topic  will  be  discussed  in  more  detail   in  Chapter  II. 
The  signal  processing  portion  of  the  relative  phase 
tern  utilizes  two  receivers,  one  for  each  antenna  pair. 

IP  outputs  of  these  receivers  are  feed  to  a  sum  and 
ference  unit  which  accurately  adds  and  subtracts  the 
rials  and  rotates  the  phase  of  either  the  sum  or  the 
ference  by  tt/2.   The  resultant  sum/difference  signals 
n  drive  orthogonal  deflection  plates  of  a  CRT  display 
t,  thus  achieving  the  effect  of  a  goniometer  rotating 
the  intermediate  frequency  of  the  receivers.   The  display 
tern  produced  by  this  technique  is  an  ellipse  with  the 
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azimuthal  line  of  bearing  along  its  major  axis.   This  signal 
processing  scheme  produces  cancellation  of  the  effects  due 
to  elevation  angle  of  arrival. 

A  system  of  this  type  is  difficult  to  build  and 
maintain  because  of  the  extreme  requirements  on  the  perfor- 
mance of  the  receivers.   Both  channels  must  be  phase  and 
gain  matched  over  the  operating  range  in  order  to  maintain 
the  amplitude  and  phase  relationship  between  the  received 
signals.   In  addition,  the  frequency  range  that  can  be 
covered  with  a  given  antenna  spacing  is  limited.   The 
receiving,  amplifying,  and  phase  comparison  portions  of  the 
system  can  be  designed  to  cover  a  wide  frequency  range. 
But  the  antenna  system  cannot  in  practice  be  made  to  cover 
a  range  wider  than  about  four  to  one  in  frequency  if 
adequate  resolving  power  without  too  many  ambiguities  is 
to  be  obtained  [3]. 

3.   Time  of  Arrival  Methods 

The  time  of  arrival  (TOA)  technique  requires  the 
measurement  of  a  given  modulation  event  at  two  known  antenna 
locations.   This  modulation  event  may  be  the  crest,  of  a 
sinusoidal  modulation  signal,  or  a  specific  point  (leading 
or  trailing  edge)  of  a  pulse  signal.   This  technique  is 
frequency  independent  and,  thus,  does  not  experience  the 
wavelength/antenna  separation  ambiguities  inherent  to  the 
relative  phase  technique.   However,  elevation  angle  compo- 
nents produce  ambiguities  just  as  in  the  relative  phase 
technique. 
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The  TOA  method  is  analogous  to  the  relative  phase 
method  in  that  antenna  baselines  are  established,  and  gain/- 
phase  matched  receivers  must  be  used  in  the  signal  processing 
section. 

k .   Doppler  Methods 

The  Doppler  principle  of  wave  motion  can  be  applied 
to  a  specially  constructed  antenna  system  to  determine  the 
angle  of  arrival  of  a  signal  [4],  If  an  omnidirectional 
antenna  is  rotated  at  a  fixed  rate  in  a  circular  path,  the 
signal  at  the  antenna  output  will  be  a  frequency  modulated 
representation  of  the  incident  signal.  The  peak  frequency 
deviation  of  this  signal  will  be: 


Af  -   ^L 
c 


where 

C  is  the  speed  of  propagation  (meters/sec), 

R  is  the  path  radius  (meters),  and 

to  is  the  rotation  rate  (radians/sec.) 

The  maximum  frequency  deviations  will  occur  when  the  antenna 
is  moving  directly  away  (negative  Af)  or  directly  toward 
(positive  Af)  the  emitter.   Although  this  property  is  not 
affected  by  elevation  angle  components  of  the  incident 
wavefront,  the  peak  deviation  frequency  will  be  diminished 
by  the  cosine  of  the  elevation  angle  and  may  reduce  the 
resolution  capability  of  the  system. 
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In  practice,  the  antenna  would  not  be  physically 
rotated;  instead  a  ring  of  omnidirectional  antennas  would 
be  constructed  and  a  capacitively  coupled  goniometer  would 
sample  each  element.   The  main  difference  between  this 
technique  and  the  Wullenweber  system  is  that  the  Doppler 
technique  samples  individual  low  gain  elements  while  the 
Wullenweber  forms  the  equivalent  high  gain  beam  pattern 
of  a  ^0-element  broadside  array. 

The  primary  disadvantage  of  this  system  is  that  it 
requires  a  circular  array  several  wavelengths  in  diameter, 
thus  precluding  its  use  in  the  lower  end  of  the  HP  frequency 
range. 
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II.   A  NEW  RDF  SYSTEM 

A.   BASIC  DESCRIPTION 

A  new  direction  finding  technique  has  been  proposed 
which  utilizes  the  finite  difference  between  arrival  times 
of  distinct  modulation  events  as  an  incident  RF  plane  wave- 
front  propagates  along  a  linear  antenna  array.   Although 
the  time  of  arrival  concept  is  not  new,  the  technique  by 
which  signals  are  processed,  coded,  and  interpreted  to 
extract  angle  of  arrival  information  is  considered  to  be 
unique.   [1] 

Basically,  the  technique  involves  the  binary  phase 
coding  and  correlation  of  the  time  separated  antenna  outputs. 
The  correlation  of  these  coded  signals  in  matched  filter 
devices  produces  a  time-compression  of  the  signals,  causing 
them  to  arrive  at  the  filter  output  at  the  same  point  in 
time.   When  this  occurs,  a  unique  correlation  output  is 
produced.   Depending  on  the  choice  of  phase  codes  utilized, 
this  unique  output  may  be  either  an  amplified  replica  of  the 
original  input  signal  (autocorrelation),  or  it  may  be  a 
distinct  null  (crosscorrelation) . 

The  characteristic  of  this  process  which  enables  the 
determination  of  arrival  angle  is  the  one-to-one  relation- 
ship between  AOA  and  the  particular  matched  filter  configura- 
tion which  produces  this  distinquishable  output.   In  practice, 
the  outputs  form  a  bank  of  matched  filters,  constructed  to 
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match  a  range  of  arrival  angles,  would  be  scanned  for  this 
unique  output.   The  filter  producing  that  output  would 
correspond  directly  to  a  particular  arrival  angle. 

The  technique  outlined  above  makes  use  of  a  special 
class  of  binary  codes  called  complementary  sequences.   The 
time  correlation  of  these  sequences  in  matched  filters  will 
produce  a  unique  output  which  is  free  of  the  coding  noise 
prevalent  with  all  other  discrete  binary  codes.   Before 
describing  in  detail  the  application  of  these  codes  to  the 
RDF  system,  however,  a  brief  discussion  of  general  phase 
coding  techniques  is  considered  necessary  for  continuity. 
Therefore,  these  techniques,  and  the  specific  advantages  of 
complementary  sequences,  are  presented  in  the  following 
section. 

B.   PHASE  CODING  TECHNIQUES 

Phase  coding  techniques  were  first  used  (in  communica- 
tions systems)  to  improve  the  range  capability  of  pulse 
radars.   While  not  relevant  to  the  topic  of  this  paper, 
this  utilization  of  phase  coding  provides  a  graphic  descrip- 
tion of  general  coding  techniques  and  an  example  of  their 
application  to  an  actual  system. 

Two  parameters  of  a  transmitted  radar  signal  may  be 
varied  in  order  to  increase  the  total  RP  energy  illuminating 
a  radar  target.   These  are  the  peak  pulse  power,  and  the 
pulse  duration.   The  peak  power  of  the  transmitted  pulses 
may  be  increased  by  utilizing  higher  power  amplifiers  in 
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the  transmitter  output  stage.   However,  this  has  proven 
inefficient  from  a  power  standpoint,  since  many  radar 
transmitter  tubes  are  peak  power  limited  while  operating  at 
only  a  small  fraction  of  their  average  power  capabilities. 
The  average  power  output  of  these  tubes  may  be  increased  by 
widening  the  transmitted  pulses.   However,  this  too  has  a 
drawback  because  the  range  resolution  capability  of  a  pulse 
radar  is  inversely  proportional  to  pulse  width. 

As  a  result  of  this  dilemma,  intrapulse  modulation 
techniques  were  developed  so  that  pulse  widths  could  be 
increased  without  sacrificing  range  resolution.   The  most 
common  example  of  intrapulse  modulation  is  the  "chirping" 
or  pulse  compression  system  in  which  the  RF  carrier  frequency 
of  each  pulse  is  swept  between  fixed  limits  during  its  trans- 
mission.  The  received  pulses  are  then  passed  through  a 
matched  filter  whose  response  characteristic  (delay  versus 
frequency)  is  the  exact  inverse  of  the  transmitted  pulse 
characteristic  (frequency  versus  time).   This  allows  all  the 
frequency  components  of  the  wide  reflected  pulse  to  arrive 
at  the  output  simultaneously,  resulting  in  a  much  narrower, 
amplified  pulse  signal  to  the  detector  stage.   Phase  coding 
of  the  transmitted  pulse  is  another  technique  of  intra- 
pulse modulation. 

The  phase  coding  technique  utilizes  a  discrete  process 
to  encode  pulses  so  that  each  return  pulse,  when  input  to  a 
matched  filter,  causes  discrete,  coded  segments  of  the  pulse 
to  arrive  at  the  output  simultaneously.   Each  radar  pulse  is 
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divided  into  "n"  equal  time  segments,  and  a  coded  sequence 
of  "phase  states"  is  assigned  to  these  segments.   The  states 
determine  the  phase  relationship  between  the  RF  carrier 
during  that  segment,  and  a  fixed  phase  carrier  reference. 
If  the  only  allowable  phase  states  in  each  segment  are  0°  or 
180°,  then  the  coding  is  termed  binary  phase  coding. 
Although  multiphase  coding  schemes  have  been  devised  and 
used  in  certain  applications,  only  binary  phase  coding  will 
be  discussed  here. 

A  representation  of  a  binary  phase  coded  pulse  is  shown 
in  Figure  2-1 (a).   The  RF  carrier  is  divided  into  seven 
equal  time  segments,  with  each  segment  assigned  a  plus  (+) 
or  minus  (-)  phase  state.   The  phase  state  values  are  here 
defined  as  no  carrier  phase  reversal  (+)  and  a  180°  carrier 
phase  reversal  (-).   A  unit  amplitude  is  assumed  for  the 
pulse . 

The  filter  which  provides  an  exact  match  for  this  coded 
pulse  is  shown  in  Figure  2-1 (b),  where  the  discrete  delay 
units  are  equal  to  the  pulse  segment  durations  and  the  (+) 
and  (-)  phase  reversals  occur  in  reverse  order  from  the 
transmitted  sequence  phase  states.   In  a  physically  realizable 
filter,  the  discrete  time  delays  would  be  produced  by  delay 
lines  or  electronic  delay  circuits,  and  the  phase  control 
would  be  provided  by  buffer  amplifiers  with  inverting  or 
non-inverting  outputs.   For  this  example,  the  buffer/phase 
control  devices  are  assumed  to  have  unity  voltage  gains. 
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Figure  2-2  is  a  time  diagram  showing  the  output  waveforms 
at  the  intermediate  points  "a"  through  "g",  and  the  final 
summation  output  waveform. 

There  are  two  characteristics  to  be  noted  about  this 
output.   First,  it  contains  a  narrow  main  pulse  that  is 
delayed  n-1  time  units  from  the  leading  edge  of  the  received 
pulse.   Since  this  narrow,  main  pulse  would  be  utilized  by 
the  radar  receiver's  detector  and  timing  circuits  to  deter- 
mine target  range,  this  delay  must  be  compensated  for. 
Second,  there  are  undesired  secondary  pulses  in  the  output 
waveform  which  might  appear  to  the  radar  receiver  as  excess 
noise  or  false  targets. 

There  are  other  codes,  called  pseudo-random  noise  codes, 
which  will  produce  smaller  secondary  pulses  than  those 
produced  by  the  (randomly  selected)  code  utilized  in  this 
example.   However,  all  single-sequence  coding  schemes  of 
this  type  produce  a  certain  amount  of  "coding  noise",  as 
these  secondary  pulses  are  called,  and  therefore  exhibit  a 
finite  processing  signal-to-noise  ratio. 

There  is  a  method,  however,  which  will  completely 
eliminate  the  coding  noise,  although  it  requires  the  simul- 
taneous transmission  of  two  distinct  coded  sequences.   These 
codes,  called  complementary  sequences,  were  first  utilized 
by  M.  J.  E.  Golay  in  the  field  of  multislit  spectrometry. 
Since  then,  they  have  been  utilized  in  various  radar  and 
sonar  system  designs.   Golay  mentioned  in  1961  that  these 
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sequences  might  be  applied  to  data  transmission  [5],  and 
there  has  been  some  work  toward  that  end  [6]. 

A  great  deal  of  work  has  gone  into  the  investigation  of 
properties  of  complementary  sequences,  and  the  methods  by 
which  they  may  be  generated.   However,  only  the  basic 
characteristics  of  these  sequences  will  be  presented  here. 
For  more  detailed  studies,  the  reader  is  referred  to 
References  5  and  7  listed  at  the  end  of  this  paper. 

A  complementary  sequence  pair  is  defined  to  be  a  pair 
of  equally  long  binary  sequences  which  have  the  property 
that  the  number  of  pairs  of  like  digits  (or  elements)  with 
any  given  separation  in  either  sequence  is  equal  to  the 
number  of  pairs  of  unlike  elements  with  the  same  separation 
in  the  other  sequence. 

For  example,  considering  the  following  sequences 

A:    ___  +  __  +  _ 
B:    ___+++_+ 

the  number  of  like  elements  in  A  with  separation  distances 
of  one  element  is  three: 

A: +  __  +  _ 

lluuluu    3   like   elements 

The  number  of  unlike  elements  in  B  with  separation  distances 
of  one  element  is  also  three: 

B:    ___+++_+ 

llulluu    3   unlike   elements 
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Similarly,  the  number  of  like  elements  in  A  and  the  number 
of  unlike  elements  in  B  for  all  separation  distances  are 
listed  below: 

Separation  Distance     #   Likes  in  A    #   Unlikes  in  B 

2  3  3 

3  4  1 
k                                             2  2 

5  2  2 

6  11 

7  11 

Golay  has  shown  [5]  that  the  sequence  lengths  must 
always  be  an  even  number,  and  must  be  the  sum  of  two  squares. 
For  example,  the  possible  sequence  lengths  up  to  50  are: 
2,  3,  8,  10,  16,  18,  20,  26,  34,  36,  40,  and  50.   He  has 
also  shown  [51  that  2   complementary  sequence  pairs  (some 
identical)  can  be  generated  by  performing  the  following  basic 
operations  on  a  given  sequence  pair: 

1.  Interchanging  the  sequences 

2.  Reversing  the  order  of  the  first  sequence 

3.  Reversing  the  order  of  the  second  sequence 

4.  Complementing  the  elements  of  the  first  sequence 

5.  Complementing  the  elements  of  the  second  sequence 

6.  Complementing  the  even-order  elements  of  both  sequences 

An  example  of  a  pair  of  complementary  sequence  coded 
pulses  being  processed  by  a  matched  filter  pair  is  shown 
in  Figures  2-3  and  2-4.   The  input  pulses,  utilizing  comple- 
mentary sequence  codes  of  length  eight,  are  shown  in  Figure 
2-3(a).   The  pulses  and  their  matched  filters,  shown  in 
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Figure  2-3(b),  are  labeled  'A'  and  !B'  respectively.   The 
outputs  from  these  matched  filters,  and  the  final  combined 
output,  are  shown  in  Figure  2-k .      This  final  output  waveform 
shows  the  signal  processing  advantage  of  this  coding  scheme. 
The  total  energy  content  of  both  received  pulses  has  been 
compressed  and  combined  into  a  single  narrow  pulse  with  an 
amplitude  gain  factor  of  twice  the  code  length.   In  addition, 
the  complete  cancellation  of  all  undesired  coding  noise 
indicates  that  an  infinite  processing  signal-to-noise  ratio 
can  be  obtained  with  this  technique. 

In  radar  system  applications,  this  property  allows  trans- 
mitters to  be  operated  at  high  average  power  levels  while 
utilizing  pulses  with  high  frequency  content.   Thus,  both 
transmitted  energy  and  range  resolution  goals  may  be  met. 
The  price  that  must  be  paid  for  this  advantage  is  the 
increased  complexity  of  the  signal  generation,  transmission, 
and  detection  components. 

The  signal  coding/decoding  process  illustrated  above  may 
also  be  described  in  mathematical  terms.   If  we  let  a.  and 
b.  be  the  elements  of  any  pair  of  n-element  binary  sequences 
(taking  on  values  +1  or  -1),  and  let  a.  and  b.  be  the  phase 
reversal  portions  of  their  corresponding  filters,  then  the 
output  sequences  from  the  filters  can  be  represented  by: 


cj  =  I     aiai+j 


d  .  =   E   b  .  b  .  .  . 
J    i=1   i  i+J 
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where  the  subscript  "j"  takes  on  integer  values  from  -n+1 
to  n-1  and,  by  definition: 


C-J  =  CJ 


d  .  =  d. 

-0     J 


The  summation  equations  represent  a  discrete  convolution 
of  the  phase  codes  and  their  respective  filter  codes, 
resulting  in  correlation  of  the  two  sequences.   In  the  case 
where  the  filters  have  "matched  filter"  characteristics 
(i.e.,  a.  =  a. )  then  the  sequences  are  said  to  be  auto- 
correlated.   These  statements  apply  to  any  pair  of  equal- 
length  binary  sequences.   Then  if  the  input  sequences 
utilized  are  a  complementary  coded  pair,  and  matched  filters 
are  used,  the  following  conditions  apply: 


s.  =  c.  +  d.  =  0        for  j  ?   0 
J    J    J 


s   =  c   +  d^  =  2n       for  j  =  0 
o     o     o  ° 


This  is  a  mathematical  description  of  the  final  output 

sequence  obtained  in  the  process  illustrated  above.   The 

final  summation  terms,  s.,  are  zero  for  all  values  except 

j 

j  =  0,  which  corresponds  to  the  main  pulse  output  during  the 
eighth  time  segment. 

A  second  useful  property  of  complementary  sequences  is 
that,  for  each  sequence  pair,  there  is  an  "orthogonal"  pair 
which  produces  a  unique  crosscorrelation  output.   That  is, 
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the  output  sequences  c.  and  d.  obtained  by  correlating  a 

J      J 

complementary  sequence  pair  with  its  corresponding  orthogonal 
pair  will  be  described  by  the  relation: 

s.  =  c.  +  d.  =0    for  all  values  of  j. 
J    J    J 

This  process  would  not  be  useful  for  pulse  radar  or  data 
transmission  systems,  since  the  processed  signal  has  zero 
amplitude  for  all  values  of  time.   However,  the  cross- 
correlation  output  is  applicable  to  the  RDF  system,  since  it 
provides  a  unique  correspondence  between  the  angle  of  arrival 
and  the  signal  processor  output. 

The  generation  of  an  orthogonal  sequence  from  a  given 
complementary  sequence  utilizes  four  of  the  six  basic 
operations  described  earlier  in  this  section.   The  generation 
algorithm,  shown  below  for  an  eight-element  complementary 
sequence  pair,  is  given  by  French  [6]  in  his  Master's  thesis. 

(1)  Given  a  pair  of  complementary  sequences: 

A  =  ---  +  --  +  - 
B  =  ---  +  +  +  -  + 

(2)  Perform  operation  #1  (interchange  sequences): 

A1  =___+++_+ 
B«  =_--+__+- 

(3)  Perform  operations  #2  and  #3  (reverse  the  order  of 
the  sequences) : 

A"  =  +  -■■+  +  +  --  - 
B"  =  -  + + 
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(*])   Perform  operation  #k    (complement  elements  of  A"): 
A"1  =  -  +  ---  +  +  + 

B'"  =  -  + + 

The  resultant  sequence  pair  A''',  B' ' '  is  orthogonal  to  the 
original  sequence  pair  A,B. 

An  example  of  the  crosscorrelation  process,  utilizing 
orthogonal  complementary  sequence  pairs,  is  shown  in  Figures 
2-5  and  2-6.   The  phase  coded  pulses  in  Figure  2-5 (a) 
utilize  the  eight  element  complementary  sequences  A  and  B. 
The  correlation  filters,  shown  in  Figure  2-5 (b),  utilize 
the  orthogonal  sequences  '''A'1'  and  ' ' 'BT ' T .   The  outputs 
from  these  filters,  as  shown  in  Figure  2-6,  add  to  zero  for 
all  values  of  time. 

With  this  introduction  to  complementary  sequence  coding 
and  an  application  to  signal  processing,  the  basis  is 
established  for  its  application  to  the  proposed  RDF  system. 

C.   APPLICATION  OF  COMPLEMENTARY  SEQUENCE  CODING  TO  AN 
IDEALIZED  RDF  SYSTEM 

As  previously  mentioned,  the  phase  coding  method  of 
determining  the  angle  of  arrival  depends  upon  the  time 
difference  between  antenna  responses  in  a  linear  array. 
This  time  difference  is  a  unique  function  of  the  inter- 
antenna  spacing,  and  the  arrival  angle  of  the  wavefront. 
This  is  illustrated  in  Figure  2-7(a),  in  which  a  single 
rectangular  RF  pulse  is  assumed  incident  to  a  four  element 
array . 
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The  rectangular  RF  voltage  pulses  induced  at  the  antenna 
outputs  are  separated  in  the  time  domain  by  a  shift  (At) 
given  by  the  equation: 


At  =  —  cos 
c 


where, 

d  is  the  antenna  spacing, 

c  is  the  propagation  velocity  of  the  wavefront, 
and       6  is  the  arrival  angle  of  the  Poynting  vector 

associated  with  the  wavefront. 
Figure  2-7 (b)  illustrates  the  relative  time  spacing  of  the 
induced  RF  voltages,  with  the  output  of  antenna  number  one 
chosen  as  a  reference.   For  a  signal  incident  upon  the 
antenna  system  at  an  arrival  angle  of  zero  (along  the  array 
baseline)  the  time  difference  will  be  equal  to  the  maximum 
propagation  time  between  antennas.   The  time  difference 
will  reduce  to  zero  as  the  arrival  angle  increases  to  the 
point  where  it  is  normal  to  the  array  baseline. 

The  type  of  antenna  used  in  the  array  is  not  a  critical 
factor  in  the  RDF  operation  of  the  system;  any  omnidirectional 
antenna  can  be  used.   However,  the  antennas  should  be  identi- 
cal in  their  electrical  characteristics  (pattern,  gain,  and 
impedance),  and  they  should  be  sufficiently  separated  to 
minimize  mutual  coupling. 

The  first  step  in  the  phase  coding  technique  is  the 
transformation  of  these  time  shifted  signals  into  a  pair 
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of  complementary  sequences  that  can  be  interpreted  by 
matched  filter  decoders.   The  basic  hardware  necessary  for 
this  process  is  shown  in  Figure  2-8.   Each  antenna  feeds 
a  dual  output  multicoupler  device  which  can  provide  both  a 
normal  phase  output  and  an  inverted  phase  output  (relative 
to  the  input  signal).   The  outputs  from  these  devices  are 
configured  so  that  their  phase  states  form  a  complementary 
sequence  pair.   For  this  reason,  the  array  length  must 
correspond  to  one  of  the  possible  complementary  sequence 
code  lengths  (2,  4,  8,  16,  18,  20,  etc.).   The  four  element 
array  was  chosen  here  for  simplicity. 

After  phase  coding,  the  output  signals  corresponding 
to  each  of  the  complementary  sequences  are  combined  in  two 
signal  channels  and  applied  to  a  code  interpreter.   The 
code  interpreter  consists  of  a  pair  of  correlation  filters 
which,  when  perfectly  matched  to  these  input  signals, 
produce  unique  autocorrelation  (or  crosscorrelation)  outputs. 

To  illustrate  the  phase  coding  and  decoding  processes 
in  more  detail,  the  four  element  array  is  used  with  an 
incident  wavefront  at  an  arrival  angle  of  zero  degrees. 
The  complementary  sequences  are  identical  to  those  in 
Figure  2-8,  and  an  idealized  rectangular  RF  pulse  is  assumed, 
with  a  pulse  width  equal  to  twice  the  inter-antenna  propa- 
gation  delay  (2—).   To  simplify  the  present  example,  it  is 
also  assumed  that  the  phase  reversal  and  summing  operations 
are  performed  at  the  carrier  frequency  of  the  incoming 
pulse  signal,  and  that  no  other  signals  are  present.   The 
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practical  requirements  for  signal  preselection,  and  inter- 
signal  discrimination  circuits  are  considered  in  a  later 
section  of  this  chapter. 

Based  on  the  above  conditions,  the  signals  at  the 
antenna  outputs  would  be  as  shown  in  Figure  2-9.   The  output 
of  antenna  number  one  is  again  chosen  as  a  reference,  and 
the  A  &  B  outputs  of  each  antenna  correspond  to  the  channels 
into  which  the  phase-coded  outputs  are  combined.   For 
convenience,  the  peak  carrier  amplitudes  of  the  wave  forms 
are  normalized  to  unity,  and  upward  or  downward  pointing 
arrows  are  used  to  indicate  their  respective  phase  states. 
An  upward  pointing  arrow  indicates  a  non-inverted  output, 
and  a  downward  pointing  arrow  indicates  an  inverted  output. 
When  these  individual  antenna  outputs  are  grouped  and  then 
added  in  an  ideal  linear  (and  lossless)  summer,  the  resultant 
output  waveforms  are  assumed  to  be  as  shown  in  Figure  2-10. 

It  should  be  noted  that  the  amplitudes  and  phases  of 
these  idealized  waveforms  are  based  on  the  assumption  that 
the  antenna  outputs  are  combined  with  perfect  carrier  phase 
cancellation  (180°  phase  difference)  or  perfect  phase  rein- 
forcement (0°  phase  difference).   To  assume  that  perfect 
phase  matching  occurs  for  all  possible  values  of  time  differ- 
ence (and,  therefore,  for  all  possible  values  of  relative 
phase  difference)  is  a  questionable  assumption.   That  is, 
perfect  180°  or  0°  matching  between  phase  coded  antenna 
outputs  will  occur  only  for  those  arrival  angles  which 
produce  phase  differences  that  are  multiples  of  2tt. 
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However,  the  waveforms  in  Figure  2-10  do  provide  a 
simplified  illustration  of  the  basic  phase  coding/decoding 
technique  applied  to  the  RDF  system.   In  addition,  it  will 
be  shown  in  the  following  section  that  the  relative  carrier 
phase  differences  between  individual  antenna  outputs  do 
not  affect  the  validity  of  the  overall  process  so  long  as 
the  signals  are  combined  in  a  linear  device  and  the  principle 
of  superposition  can  be  applied. 

Returning  then  to  this  simplified  example,  the  coded 
and  summed  signals  occupying  channel  A  and  channel  B  are 
each  applied  to  correlation  filters  which  are  configured 
as  shown  in  Figure  2-11.   The  filters  consist  of  discrete 
delay  devices  (or  tapped  delay  lines)  and  dual  output 
buffer  devices  with  inverted  and  non-inverted  outputs.   The 
dual  outputs  of  the  buffer  units  are  configured  so  that  both 
autocorrelation  and  crosscorrelation  of  the  coded  input 
signals  can  be  performed.   For  perfect  filter-to-signal 
matching,  the  filter  time  delays  must  equal  the  inter- 
antenna  propagation  time  (—  cos  8).   For  this  example, 
(6  =  0°)  perfect  matching  requires  that  the  time  delays 
equal  — .   The  autocorrelation  and  crosscorrelation  filter 
outputs  are  shown  in  Figure  2-12  for  the  perfect  matching 
case.   As  in  Figure  2-9,  upward  and  downward  pointing  arrows 
indicate  the  signal  phase  states,  and  the  amplitudes  corres- 
pond to  peak  values.   When  the  corresponding  (autocorrela- 
tion and  crosscorrelation)  outputs  of  channels  A  and  B  are 
combined  in  a  linear  device,  the  final  output  waveforms  will 
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be  as  shown  in  Figure  2-13.   The  autocorrelation  output, 
for  perfect  filter-to-signal  matching,  is  an  exact  replica 
of  the  original  RF  pulse,  with  a  processing  gain  equal  to 
twice  the  code  length.   The  crosscorrelation  output  is 
everywhere  zero  because  of  the  complete  phase  cancellation 
of  all  signal  components. 

A  computer  model  of  this  process  has  been  used  [1]  to 
determine  the  effect  of  correlation  filter  mismatch  upon 
the  autocorrelation  and  crosscorrelation  output  waveforms. 
An  important  result  of  that  study  was  the  determination 
that  the  crosscorrelation  output  is  uniquely  zero  for  the 
perfectly  matched  case,  and  that  a  very  slight  amount  of 
mismatch  will  cause  a  rapid  increase  in  the  crosscorrelation 
output.   It  was  concluded,  then,  that  the  distinct  null 
effect  produced  by  the  crosscorrelation  process  could  be 
utilized  to  accurately  determine  the  angle  of  arrival  of 
incident  signals. 

The  computer  model  utilized  in  that  study  employed  the 
idealized  signal  model  that  was  assumed  above,  and  it  also 
assumed  a  noise-free  signal  environment.   In  order  to  more 
fully  evaluate  the  performance  of  the  complementary  sequence 
phase  coding  process  for  this  paper,  the  computer  model  was 
modified  to  investigate  the  effects  of  background  noise  and 
various  signal  parameters.   In  addition,  the  model  was 
utilized  to  investigate  the  effects  of  array  length  upon  the 
overall  process.   A  description  of  the  model,  and  the  results 
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and  conclusions  reached,  are  presented  in  Chapter  III  of 
this  paper. 

D.   SINUSOIDAL  ANALYSIS 

The  following  sinusoidal  analysis  of  the  coding/decoding 
process  will  show  that  the  basic  technique  schematically 
described  in  the  previous  section  is  valid  for  all  arrival 
angles  (between  0°  and  90°),  and  that  a  0°  or  180°  phase 
relationship  between  the  antenna  outputs  is  not  a  require- 
ment for  proper  operation  of  the  system. 

First  of  all,  it  is  assumed  that  the  physical  constraints 
on  the  system  (array  size,  phase  coding  sequences,  matched 
filter  configuration,  and  signal  parameters)  are  identical 
to  those  specified  in  the  previous  example.   However,  for 
the  purpose  of  this  analysis,  the  arrival  angle  may  take  on 
any  arbitrary  value  between  0°  and  90°. 

Therefore,  the  time  differences  between  antenna  responses 

will  be 

At  =  —  cos  8 
c 

and  the  relative  phase  differences  will  be 

A$  =  2tt  y  cos  6  =  oj  —  cos  0   , 

A  C 

where  w  is  the  carrier  frequency  of  the  incident  signal. 
Now,  if  antenna  number  one  is  chosen  as  a  time  and  phase 
reference,  then  the  signal  voltages  induced  in  the  four 
antennas  may  be  represented  by  the  equations: 
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v-,(t)  =  [u(t  )-u(t-2— )  ]cos  wt  =  u,  (t)cos  cat 


v2(t)  =  [u(t-At)-u(t-2--At)]cos(a)t-A<j>)  =  u2(t  )cos  (u)t-A<J>) 


VQ(t)  =  [u(t-2At)-u(t-2--  2At)]cos(u)t-2A<{>)  =  uq(t  )cos(a)t-2A<|>) 
j  c  j 


Vh(t)  =  [u(t-3At)-u(t-2-- 3At)]cos(u>t-3A<(>)  =  Uj, (t )cos (wt-3A(J>) 


These  equations  represent  the  signal  waveforms  shown  in 
Figure  2-7.   The  gating  functions  u, (t)  through  u^(t)  define 
the  leading  and  trailing  edges  of  the  signal  pulses,  which 
are  of  duration  2  —  .   These  signals  are  then  phase  coded 
using  complementary  sequences  and  combined  in  a  pair  of 
linear  devices  so  that  the  resultant  outputs  are  represented 
by: 


Channel  A:   v«  =  +v,  +  v~  +  v,  -  v^ 


Channel  B:   vR  =  +v,  +  v„  -  v-  +  vn 


The  composite  signals  v.  and  vR  are  then  input  to  auto- 
correlation and  crosscorrelation  filters  with  discrete  time 
delay  segments  equal  to  —  cos  0.   Since  the  time  delays  are 
equivalent  to  negative  phase  shifts  (-A<f>)  equal  to  -oj(—  cos  6), 
then  the  four  correlation  outputs  can  be  represented  by  the 
series : 


40 


Channel   A    (autocorrelated) : 

VA   =    "VA   +   va/"A())      +   vA/"2A(f)      +   vA/"3Al 


Channel  B  (autocorrelated): 

VB  =  +VB  "  Vp/"A(|)   +  vd/-2A(^)   +  vE/-3A<{) 


Channel  A  (crosscorrelated) : 

VA  =  "VA  "  VA7"A(1)   +  VA/"2A(')   "  va/"3A(|) 


Channel  B  (crosscorrelated): 

VB  =  +vb  +  Vp/-A^)   +  vD/-2A(ft   "  vc/-3A^ 


where  vy/-nA(f)  represents  the  composite  voltage  (v.  or  vR) 
phase  shifted  by  an  angle  -nA<j>.   For  example,  expanding 
v./-nA<t>  produces  the  following  equation: 

v  ./-nAcj)  =  v..  /-nA<j)   +  v0/-nAcf)   +  v0/-nA<{)   -  V|./-nA(j> 


A  general  expression  for  the  series  terms  is 


vk/-nA(t>  =  {uk(t)  cosCwt  -  (k-l)A<J>]J  /-nAcj)    Cn-0*1*2*3) 


where  u,  (t)  =  u[t  -  (k-l)At]  -  u[t  -  2  --  (k-l)At]  . 

rC  C 


Since  the  phase  delay  /-nA(fr   is  equivalent  to  the  time  delay 
nAt,  the  gating  function  u.(t)  is  delayed  by  nAt,  and  becomes 


ufc(t+nAt)  =  u[t  -  (k+n-l)At]  -  u[t  -  2 --  (k+n-1)  At  ]  =  uk+n(t) 


kl 


Therefore,  the  general  series  term  may  be  expressed  as: 
v,  /-n A<f>   =  u,   (t)  cos[u)t  -  (k+n-l)A<(>] 

Now,  applying  these  terms  to  the  autocorrelation  and  cross- 
correlation  series: 

(k=l)  (k=2)  (k=3)  (k-4) 

V  = -u,  cos  (cot )  -Upcos(oot-A4))    -u.-.cos(cot-2A4>)  +  Uj.cos(tot-3A<|>)       (n=0) 

+UpCOs  (a)t-A<J> )  +  u-cos  (wt-2A<J> )  +  Uj, cos  (cot-3Acj> )  -  u^cos  (cot-4A<}> )  (n=l ) 
+u~cos  (cot-2A({) )  +  Uu cos  (wt— 3A<t> )  +  Uj-cos  (cot-4 A(J> )  -  u,-cos (cot-5A<f> )  (n=2 ) 
-Uj,cos(a)t-3A<|>)  +Ut-cos(ajt-^A(|))  +  u^cos  (u>t-5A<|> )  -u7cos(cot-6A<}))       (n=3) 

= -u-,  cos  (cot )  +  iucos  (o)t-2A<J> )  +4ui,cos(cot-3A<{))  +uI-cos(a)t-ilA(J>)  -UyCOs(wt-6A<J)) 

(k=l)  (k-2)  (k=3)  (k-4) 

V  =+u-,cos(cot)  +UpCos(cot-A<}))    -u~cos(wt-2A<|>)  +  U|,cos(tot-3A<t>)        (n=0) 

-UpCOs(cot-A<fO  -u^cos(cot-2A4i) +ui,cos(tot-3A(})) -uc-cos(cot-iiAc{))  (n=l) 
+iucos(cot-2A<J>)  +UhCos(oot-3A(J))  -uc-cos(cot-^A({))  +  u,-cos(cot-5A4))  (n=2) 
+U|.cos(cot-3A<J>)  +  uc-cos(cot-iJAct>)  -u,-cos(cot-5A4>)  +u„cos(cot-6A({))       (n=3) 

=   u..  cos  (cot)  -iucos(cot-2A<f>)  +  4uhCOs(cot-3A<f>)  -ut-cos(cot-ilA4))  +u7cos(tot-6A<})) 

(k=l)  (k-2)  (k=3)  (k-4) 

t 

V  =-uTcos(cot)  -UpCOs(wt-A<}))    -u^cos(tot-2A({))  +  u.nCos(u>t-3A<J>)       (n=0) 

-UpCos(cot-A4))  -u„cos(tot-2A4))  -UhCos(wt-3A4>)  +uc-cos(cot-iJA<|>)  (n=l) 
+u~cos(cot-2A<t>)  +  u^cos(cot-3A<t>)  +u,-cos(tot-4A({))  -  UgCos(tot-5A<}>)  (n=2) 
-U2,cos(cot-3A4))  -ut-cos(tot-ilA(|))  -UgCos(tot-5A<}>)  +u7cos(wt-6A<j>)       (n=3) 

=  -u1cos(tot)-2u2cos(a)t-A4))-UoCos(a^-2A^)+UcCOs(cot-^Ac{))-2u£COs(tot-5Ac()) 
+u7cos(tot-6A<j>). 
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(k=l)  (k=2)  (k=3)  (k=4) 

V  =+u-,cos(u)t)  +UpCos(a)t-A(f))    -u_cos(a)t-2Ac{>) +U2.cos(ajt-3A<iO  (n=0) 

+UpCos(wt-A({))    +UoCos(cjt-2A4>) -UhCOs(cot-3A(j)) +ut-Cos(ajt~iIA(j))  (n=l) 

+u~cos(wt-2A4>)  +UhCos(cot~3A<J))  -u,_cos(u)t-*IAc{>)  +  U£-cos((jot-5Acf))  (n=2) 

-U|.cos(a)t-3Acj))  -uucosCwt-^Acf))  +  u,-cos(u)t-5A<j))  -  iucos((jot-2Acj))  (n=3) 

=  u-.cos(uyt)  +  2UpCos(ayt-Ac}>)  +  UoCos(wt-2AcJ))  -  LucosCtat-^Atj)) 
+  2ivcos(wt-5Acj))  -  u7cos(wt-6Ac{)). 

Finally,  when  each  channel's  autocorrelation  and  cross- 
correlation  terms  are  combined,  the  resulting  outputs  will 
be: 

Vq  =  VA  +  VB  =  Buncos  (o)t-3A«J)) 

=    8Ju[t-3At]    -   u[t  -  2  |-3At]}cos(a)t-3A(j)) 

Vo    =   \   +   VB   =    ° 

These  equations  represent  the  final  output  waveforms  shown 
schematically  in  Figure  2-13.   The  autocorrelation  output  is 
an  exact  replica  of  the  incoming  pulse  signal,  and  is 
amplified  by  a  factor  of  twice  the  code  length.   The  gating 
function  Uj.(t)  corresponds  to  the  result  in  Figure  2-13 
in  that  the  leading  edge  of  the  output  pulse  is  delayed 
3At  time  units  from  the  initial  intercept  of  the  incident 
pulse.   Also,  as  in  the  graphical  example,  the  final 
crosscorrelation  output  is  everywhere  zero.   Therefore,  it 
is  concluded  that  the  phase  coding  process  is  valid  for  all 
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arrival  angles  between  0°  and  90° ,  that  it  provides  a  signal 
processing  gain  (in  the  autocorrelation  mode)  of  twice  the 
chosen  code  length,  and  that  it  provides  a  crosscorrelation 
output  that  is  everywhere  zero. 

E.   ENGINEERING  CONSIDERATIONS 

In  order  to  apply  the  basic  phase  coding  technique 
described  above  to  a  physically  realizable  RDF  system, 
certain  practical  problem  areas  must  be  considered.   During 
the  investigation  and  evaluation  of  this  technique,  it 
became  obvious  that  several  factors  have  a  direct  effect 
upon  the  operation  of  the  proposed  system.   Specifically, 
those  factors  that  were  considered  most  important  are: 
(1)  possible  sources  of  error  resulting  from  angular 
ambiguities,  (2)  the  requirement  for  signal  selection  and 
discrimination,  and  (3)  the  effect  of  different  modulation 
types  on  the  correlation  processes.   The  constraints  and 
physical  requirements  which  these  factors  place  upon  the 
design  of  such  a  system  will,  in  large  part,  determine  the 
feasibility  of  utilizing  this  type  of  RDF  system  over  some 
other  type . 

The  first  two  factors  (angular  ambiguity  and  signal 
preselection)  are  discussed  and  evaluated  in  this  section. 
The  actual  design  requirements  that  must  be  imposed  to 
satisfy  these  two  areas  are  straightforward,  and  are 
considered  to  be  technically  feasible.   However,  the  overall 
practicality  or  economic  feasibility  of  these  requirements 
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will  depend,  in  part,  upon  the  particular  application  of 
the  system,  and  are  not  addressed  in  this  paper.   Modulation 
effects  are  best  described  and  evaluated  by  the  computer 
model  used  to  simulate  the  phase  coding  and  decoding  portion 
of  the  system.   This  discussion  is  presented  in  the 
following  chapter. 

1.   Angular  Ambiguities 

Until  now,  the  determination  of  arrival  angle  has 
been  assumed  to  be  a  two-dimensional  problem,  in  which  one 
single-variable  equation  is  evaluated  to  yield  a  single 
valued  unambiguous  solution.   However,  in  most  environmental 
situations,  the  direction  of  arrival  of  an  incident  wave- 
front  will  generally  be  composed  of  both  azimuthal  and 
elevation  angle  of  arrival  components.   The  additional 
consideration  of  elevation  angle  introduces  a  third  spatial 
dimension  to  the  previous  two  dimensional  problem.   That  is, 
the  elevation  angle  introduces  a  second  variable  to  the 
time  delay  equation,  so  that  At  =  —  cos  G  cos  a,   where  a 
is  the  elevation  angle  of  arrival.   This  additional  variable 
precludes  the  unambiguous  determination  of  azimuthal  angle 
from  the  single  time  difference  equation. 

Furthermore,  quadrantal  ambiguities  exist  because 
of  the  periodic  nature  of  the  cosine  function.   For  example, 
azimuth  arrival  angles  of  ^5°,  135°,  225°,  and  315°  would 
produce  equal  magnitude  time  differences  between  adjacent 
antenna  responses.   Fortunately,  some  differentiation 
between  these  responses  is  possible  by  knowing  the  algebraic 
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sign  of  the  time  differences.   This  differentiation  can  be 
accomplished  by  a  system  using  the  phase  coding  technique 
because  the  chronological  order  in  which  the  individual 
antennas  respond  to  the  incident  wavefront  will  uniquely 
determine  the  composition  of  the  phase  coded  sequences. 
For  example,  two  incident  wavefronts  whose  azimuthal  arrival 
angles  differ  by  180°  will  create  distinctively  different 
phase  coded  sequences.   Correlation  filter  pairs  which  are 
matched  to  either  one  of  these  sequences  will  not  be  matched 
to  the  other. 

This  can  be  demonstrated  using  the  example  illustrated 
in  Figures  2-8  and  2-11,  with  the  rectangular  pulse  approach- 
ing from  l80°  instead  of  0°.   Assuming  the  same  phase  coding 
of  the  antenna  signals,  the  combined  antenna  outputs  for 
the  180°  AOA  signal  will  be  as  shown  in  Figure  2-l4.   When 
those  signals  are  input  to  the  autocorrelation  and  cross- 
correlation  filters  (which  remain  matched  for  0°  AOA  signals) 
the  output  waveforms  will  be  as  shown  in  Figure  2-15.   Com- 
bining the  respective  autocorrelation  and  crosscorrelation 
outputs  from  the  channel  A  and  channel  B  filters,  the  final 
output  waveforms  will  be  as  shown  in  Figure  2-l6.   These 
outputs  indicate  a  severe  mismatch  between  the  coded  signals 
and  the  correlation  filters.   Therefore,  there  is  no  sense 
ambiguity  in  the  response  of  the  system.   In  an  actual  RDF 
system  utilizing  this  technique  it  would  be  possible  to 
establish  a  bank  of  correlation  filters  which  would  respond 
(uniquely)  to  a  full  180°  range  of  azimuthal  arrival  angles. 
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However,  there  are  still  ambiguities  created  by 
azimuthal  angles  which  produce  equal  (sign  and  magnitude) 
time  differences.   These  ambiguities,  and  the  ambiguities 
created  by  elevation  angles  of  arrival  must  be  resolved  by 
an  interferometer  technique.   A  pair  of  linear  arrays, 
situated  in  the  ground  plane  and  separated  by  a  known  fixed 
angle,  will  yield  distinct,  simultaneous  time  difference 
equations  related  to  their  respective  baselines.   A  brief 
analysis  of  the  two-array  system,  which  is  represented  in 
Figure  2-17,  will  show  that  all  spatial  ambiguities  can  be 
resolved. 

The  time  difference  equations  which  describe  the 
(orthogonal)  arrays  in  Figure  2-17  are: 


(At)..  =  -  cos  Q      cos  a 
1    C      1 


(At)„  =  —  cos  60  cos  a 
2    c      £ 


where  the  subscripted  time  shifts  (At).,  and  (At)?  are  intro- 
duced by  the  corresponding  linear  arrays,  and  a  represents 
the  fixed  elevation  angle,  which  may  take  on  any  values 
between  0°  and  90°.   The  azimuth  angles  8,  and  8p  are  related 
by  the  orientation  of  the  two  linear  arrays.   For  example, 
in  the  case  of  the  orthogonal  arrays  in  Figure  2-17; 


e2  =  e  +  270°. 
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These  three  equations  will  yield  unambiguous  solutions 
for  the  azimuth  and  elevation  angles  if  the  amplitude  and 
sign  of  (At),  and  (At)?  can  be  determined.   The  uniqueness 
of  these  equations  can  be  shown  using  the  four  incident 
wavefronts  represented  in  Figure  2-17  by  the  vectors  I,  II, 
III,  and  IV.   An  arbitrary  elevation  angle  (a)  is  assumed 
for  all  four  wavefronts,  and  their  azimuth  angles  are  chosen 
as  30°,  150°,  210°,  and  330°  relative  to  the  axis  of  antenna 
number  one.   The  time  difference  equations  corresponding  to 
these  four  wavefronts  are  as  follows: 


Case  I.     (6   =  30°,  Q?   =  300°) 


(At),  =  —  cos  9_  cos  a   =  +0.87  —  cos  a 

1  c  1  c 

(At)0  =  —  cos  G„  cos  a  =  +0.50  —  cos  a 

2  c  2  c 


Case 


II.    (9  =  150°,  e2  =  60°) 


(At),  =  —  cos  0,  cos  a  =  -O.87  —  cos  a 

1  c  1  c 

(At)0  =  —  cos  0~  cos  a  =  +0.50  —  cos  a 

2  C  2  C 


Case  III:   (61  =  210°,  ©2  =  120°) 


(At),  =  —  cos  9,  cos  a  =  -O.87  —  cos  a 

J-  w  JL                                                  w 

(At)0  =  —  cos  9n  cos  a  =  -0.50  —  cos  a 

2  c  2               c 
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Case  IV:    (6   =  330°,  92  =  2^0°) 


(At),  =  —  cos  0,  cos  a  =  +0.87  —  cos  a 

1  c  1  c 

(At)~  =  —  cos  6~  cos  a  =  -0.50  —  cos  a 

2  c  2  c 


Although  the  corresponding  time  shifts  introduced 
by  each  of  the  wavefronts  have  equal  magnitudes,  their 
algebraic  signs  are  unique  pair  combinations.   Therefore, 
differentiation  between  these  (otherwise  ambiguous)  wave- 
fronts  can  be  performed  once  the  algebraic  sign  of  each  time 
shift  is  known.   That  is,  each  pair  of  time  difference 
values  corresponds  uniquely  to  a  given  azimuth  angle  and 
elevation  angle  combination. 

It  has  already  been  shown  that  the  "sense"  and 
therefore,  the  algebraic  sign  of  the  time  shifts  can  be 
determined  by  the  matched  filter  correlation  process. 
Therefore,  the  two  array  system  suggested  here  is  theoret- 
ically capable  of  determining  the  parameters  necessary  to 
locate  the  incident  wavefront  in  three  dimensional  space. 
However,  the  design  of  a  specific  mechanism  to  accomplish 
this  is  not  attempted  in  this  paper.   The  design,  modeling, 
and  evaluation  of  a  physically  realizable,  three-dimensional 
RDF  system  utilizing  phase  coded  linear  arrays  is  considered 
a  sufficient  task  for  additional  thesis  work. 
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2.   Signal  Pre-Selection 

As  in  any  radio  frequency  receiving  systems,  the 
radio  direction  finder  must  be  capable  of  electronically 
sorting  through  a  range  of  operating  frequencies  in  order 
to  select  particular  signals  of  interest.   In  most  RDF 
systems  currently  in  use,  pre-selection  and  identification 
of  signals  is  accomplished,  either  by  a  human  operator  or 
automatically,  before  the  line  of  bearing  is  obtained  for 
that  signal.   The  idealized  RDF  system  presented  in  this 
paper  has  thus  far  been  limited  to  operation  in  a  single- 
signal  environment.   However,  in  order  to  utilize  the  phase 
coding  technique  in  a  realizable  system  capable  of  operating 
in  a  dense  RF  environment,  signal  preselection  and  amplifi- 
cation functions  must  be  incorporated  into  the  signal 
processing  portion  of  the  system. 

The  pre-selection  method  which  is  used  almost 
universally  in  communications  receivers  is  the  superhetero- 
dyne technique.   That  is,  the  desired  signal  frequency  is 
converted  by  means  of  a  non-linear  mixing  process  to  some 
intermediate  frequency  (IF)  at  which  narrowband  high  gain 
amplifiers  may  be  used.   In  some  cases  the  frequency  con- 
verter circuit  is  preceded  by  one  or  more  tuned  RF  amplifiers 
to  provide  coarse  frequency  selection,  and  to  improve  the 
receiver  noise  figure.   However,  the  overall  receiver 
selectivity  is  determined  almost  entirely  by  the  IF  circuits 
employed. 
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In  order  to  apply  the  superheterodyne  principle  to 
the  proposed  RDF  system,  a  decision  must  be  made  concerning 
where  the  frequency  conversion  should  occur  in  the  signal 
coding  and  decoding  process.   This  decision  will  influence 
the  type  of  circuits  and  components  that  would  be  required 
in  the  implementation  of  such  a  system,  and  in  some  cases 
may  determine  whether  or  not  the  system  is  operationally 
feasible.   Basically,  the  conversion  must  occur  at  a  stage 
of  signal  processing  where  it  will  not  affect  the  RDF  opera- 
tion of  the  system.   There  are  several  locations  in  the 
signal  path  where  frequency  conversion  and  envelope  detection 
might  take  place. 

The  first  possibility  involves  the  phase  coding 
and  correlation  of  all  input  signals  at  the  carrier  fre- 
quency, with  signal  selection,  conversion,  and  envelope 
detection  at  the  output  of  the  matched  filters.   The  detected 
outputs  of  the  matched  filters  would  then  be  scanned  for  a 
unique  crosscorrelation  output.   This  approach  has  been 
utilized  in  the  discussion  and  analysis  of  the  phase  coded 
system  presented  thus  far  in  this  paper.   The  sinusoidal 
analysis  in  the  previous  section  has  shown  that  phase  coding 
and  matched  filter  correlation  of  signals  can  be  performed 
at  the  carrier  frequency  for  all  arrival  angles. 

The  primary  disadvantage  of  this  approach  is  that 
the  signal  processing  is  performed  at  very  low  signal  levels. 
Accurate  180°  phase  inversion  and  amplitude  balance  of 
antenna  outputs  over  a  wide  range  of  frequencies  would 
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seem  to  be  a  very  difficult  requirement,  considering  that 
signal  levels  are  typically  in  the  microvolt  range.   These 
signal  levels  would  be  further  reduced  by  the  correlation 
filters,  since  most  time  delay  devices,  such  as  coaxial 
cable,  electromechanical  delay  lines,  and  acoustic  surface 
wave  devices,  are  passive,  and  introduce  signal  losses. 

In  addition,  it  is  not  clear  what  effect  many 
different  signals  at  various  frequencies  and  arrival  angles, 
will  have  on  the  phase  coding  and  correlation  processes. 
The  possibility  of  signal  interactions  due  to  circuit 
non-linearities,  which  are  always  present  in  active  devices, 
indicates  that  ambiguous  results  or  spurious  responses  might 
be  introduced  by  this  wide-band  approach. 

This  approach  does  have  the  advantage  however,  that 
only  one  receiver  is  required.   That  is,  the  various  antenna 
signals  are  eventually  combined  into  a  single  channel  which 
could  then  be  input  to  a  standard  communications  receiver. 
In  this  case  the  scanning  process  could  be  performed  by 
electronically  sampling  the  matched  filter  outputs  and 
feeding  the  samples  to  the  single  receiver. 

A  second  possible  configuration  of  the  signal  pro- 
cessor involves  the  frequency  conversion  of  signals  at  each 
of  the  antenna  outputs.   With  this  procedure,  the  low  level 
RF  signals  are  converted  to  IP  signals  before  the  phase 
coding/decoding  process  occurs.   In  this  case,  the  single- 
signal  situation  may  again  be  assumed,  with  the  coding/decoding 
being  performed  at  the  intermediate  frequency  instead  of 

the  carrier  frequency. 
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This  approach  would  seem  to  offer  several  advantages 
over  the  previous  wide-open  front  end  approach.   First  of 
all,  signal  preselection  before  phase  coding  and  correlation 
processing  will  eliminate  the  possibility  of  interaction 
between  signals.   Second,  the  noise  figure  of  the  system 
will  be  established  by  the  individual  array  receivers,  so 
the  losses  introduced  by  the  correlation  process  will  have 
negligible  effect  upon  the  overall  sensitivity  of  the 
system.   Finally,  since  the  phase  coding  and  correlation 
process  will  be  performed  on  high  level,  fixed  frequency 
signals,  more  efficient  narrowband  circuits  can  be  employed. 

The  primary  disadvantage  of  this  system  is  that  it 
requires  the  use  of  multiple  receivers,  which  must  be  gain 
and  phase  matched  over  their  total  frequency  range.   However, 
it  will  be  shown  in  the  next  chapter,  that  the  use  of  large 
arrays  is  not  necessary  for  good  RDF  resolution  capability. 
Therefore,  the  use  of  short  (2  or  k    element)  array  with 
this  approach  might  result  in  a  satisfactory  trade-off 
between  system  performance  and  complexity.   This  approach, 
utilizing  individual  receivers  at  each  antenna,  and  phase 
coding  and  correlation  of  signals  at  a  fixed  intermediate 
frequency,  is  employed  in  the  system  model  evaluated  (by 
computer  simulation)  in  the  following  chapter. 
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Figure  2-1 ,  Phsse  Coded  Input  Pulse  and  Matched  Filter 
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Figure  2-2,  Matched  Filter  Output  Waveforms 
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Figure  2-3,  Complementary  Sequence  Phase  Coded  Pulses  and  Matched 

Filters 
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Figure  2-4,  Matched  Filter  Output  Waveforms 
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Figure  2-5,  Complementary  Sequence  Phase  Coded  Pulses  and  Matched 
Crosscorrelation  Filters 
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Figure  2-6,  Crosscorrelation  Filter  Output  Waveforms 
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Figure  2-7,  Four  Element  Antenna  Array  and  Induced  Signals 
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Figure  2-8,  Antenna  Output  Phase  Coding  Scheme 
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Figure  2-9,  Antenna  Output  Waveforms 
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Figure  2-10,  Combined  Output  Waveforms 
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Figure  2-11,  Four  Element  Correlation  Filters 
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Figure  2-12,  Correlation  Filter  Output  Waveforms 
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Figure  2—13,  Final  Autocorrelation  and  Crosscorrelation  Outputs 
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Figure  2-1 4.,  Combined  Antenna  Outputs  for  180  Arrival  Angle 
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Figure  2-15,  Correlation  Filter  Outputs  for  180°  Mismatch 
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Figure  2-16,  Final  Autocorrelation  and  Crosscorrelation  Outputs 

For  a  180  Mismatch 
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Figure  2-17,  Two  Array  Interferometer 
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III.   SYSTEM  MODELING 

In  order  to  predict  the  effectiveness  of  the  phase  coded 
RDF  technique  in  an  actual  signal  environment,  certain 
aspects  of  its  performance  should  be  evaluated.   As  in  any 
RDF  system,  the  most  important  parameters  by  which  its 
usefulness  can  be  judged  are  its  angular  resolution  capabil- 
ity, and  its  performance  in  a  noisy  environment.   It  has 
already  been  shown  that  an  idealized  two  array  system  is 
possible  which  can  eliminate  the  gross  bearing  errors  due 
to  angular  ambiguities.   However,  the  basic  performance 
parameters  which  will  govern  its  precision  are  those  inherent 
to  the  individual  linear  arrays.   Specifically,  the  bearing 
accuracy  and  resolution  capability  of  the  system  is  funda- 
mentally dependent  on  the  precision  of  the  signal  processing 
technique  used  to  determine  time  of  arrival  differences. 

Therefore  a  computer  simulation  of  the  signal  encoding 
and  decoding  portion  of  the  RDF  system  was  employed  to 
evaluate  its  performance  under  noisy  signal  conditions,  and 
to  determine  the  effects  of  varying  array  lengths  and  signal 
parameters.   The  program  was  originally  developed  by  R.  C. 
Todaro  in  his  Master's  thesis  [1]  to  test  the  uniqueness  of 
the  correlation  outputs  for  varying  degrees  of  correlation 
filter  mismatch.   The  original  signal  model  used  in  the 
simulation  was  an  idalized  rectangular  RF  pulse  of  fixed 
duration,  and  a  noise-free  environment  was  assumed.   The 
output  from  this  simulation  consisted  of  printer  plots  of 


71 


autocorrelation  and  crosscorrelation  outputs  for  various 
arrival  angles,  and  for  varying  degrees  of  correlation 
filter  mismatch. 

The  correlation  results  obtained  from  the  original  model 
indicated  that  no  inherent  ambiguities  existed  for  arrival 
angles  within  any  given  quadrant  of  the  array  ground  plane. 
In  addition,  the  crosscorrelation  plots  obtained,  using  a 
l6-element  array  model,  indicated  that  the  angular  resolution 
of  the  system  was  limited  only  by  the  precision  of  the 
matched  filter  delay  units.   That  is,  if  delay  lines  could 
be  physically  realized  so  that  their  time  delays  corre- 
sponded to  each  incremental  degree  of  arrival  angle,  then 
the  resolution  of  the  system  would  be  ±0.5°. 

A.   MODEL  DESCRIPTION 

For  the  purpose  of  this  paper,  the  program  was  modified 
to  produce  continuous  fine-grain  autocorrelation  and  cross- 
correlation  output  plots  from  a  CALCOMP  plotter  for  varying 
noise  levels,  array  lengths,  and  input  signal  parameters. 
The  system  model  is  identical  to  the  graphical  illustration 
used  in  the  previous  chapter.   A  single  linear  array  was 
used,  with  the  azimuthal  angle  of  arrival  measured  from  the 
baseline  of  the  array.   The  signal  model  was  modified  as 
shown  in  Figure  3-1,  to  investigate  the  effects  of  input 
signal-to-noise  ratio  on  the  correlation  outputs.   A 
constant  noise  level,  corresponding  to  the  average  (or  RMS) 
level  in  an  actual  signal  environment,  was  chosen. 
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Figure  3-1.   Idealized  Signal  Model 


The  signal-to-noise  ratio  was  defined  as  the  ratio  of  the 
peak  signal  voltage  to  the  RMS  noise  voltage. 

The  model  parameters  that  are  variable,  either  by  choice 
of  specific  program  data  cards,  or  by  internal  adjustments 
to  program  instructions,  are  as  f ollows : 

1.  Array  length  (2,4,8,10,16,  and  18  may  be  selected  by 
data  cards,  larger  arrays  require  re-dimensioning  of 
program  storage). 

2.  Arrival  angle  (selected  by  data  cards) 

3.  Input  signal-to-noise  ratio  (selected  by  two  program 
instructions) . 

4.  Choice  of  autocorrelation  or  crosscorrelation  outputs 
(selected  by  data  cards). 

5.  Complementary  sequences  (any  valid  complementary  code 
pair -may  be  loaded  via  data  cards,  the  filter  codes 
for  autocorrelation  or  crosscorrelation  processes  are 
automatically  derived  by  the  program) . 
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6.  Filter/Signal  matching  coefficient,  D  (the  ratio  of 
individual  correlation  filter  delay  to  maximum 
propagation  delay  At,  selected  by  data  card). 

7.  Input  pulse  width  (selected  by  program  instruction). 
A  detailed  discussion  of  the  program  operation,  format 

of  input  data,  and  an  instruction  listing  are  included  in 
the  appendix. 

B.   MODEL  OUTPUT  RESULTS  AND  CONCLUSIONS 

The  output  plots  obtained  from  this  simulation  are 
divided  into  three  specific  groups  in  order  to  show  the 
effects  of  individual  parameter  changes  on  the  overall 
operation  of  the  model.   These  groups  are  outlined  below, 
and  are  referenced  to  specific  figure  numbers,  so  that  the 
related  discussion  and  conclusions  may  be  presented  together 
at  the  end  of  this  section.   All  pertinent  parameters  are 
listed  in  the  heading  of  each  of  the  plots. 
1.   Group  1,  Figures  3-2  through  3-^1 

These  plots  demonstrate  the  effect  of  varying  the 
input  signal-to-noise  ratio  on  autocorrelation  and  cross- 
correlation  outputs.  SNR  takes  on  values  of  infinity  (no 
noise),  10:1,  5:1,  and  2:1. 

Fixed  parameters:   Arrival  Angle  =  60.0° 

Array  Length  =  16 

Pulse  Duration  =  2.0  times  maximum 

inter-antenna 
propagation  time. 

For  each  combination  of  SNR  and  fixed  parameters  the  matching 

coefficient  (D)  takes  on  values  of  0.900,  0.950,  1.000,  1.050 

and  1.100. 
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2.  Group  2,  Figures  3-^2  through  3-56 

These  plots  demonstrate  the  effect  of  varying  array 
length  on  autocorrelation  and  crosscorrelation  outputs. 
Array  lengths  of  10,  8,  and  4  elements  are  utilized. 

Fixed  parameters:   Arrival  angle   =  60.0° 

Pulse  Duration  =  2.0  times  maximum 

inter-antenna 
propagation  time. 

For  each  combination  of  array  length  and  fixed  parameters, 

the  matching  coefficient  takes  on  values  of  0.900,  0.950, 

1.00,  1.050  and  1.100. 

3.  Group  3,  Figures  3-57  through  3-80 

Effect  of  varying  pulse  width/antenna  separation 
ratio  on  autocorrelation  and  crosscorrelation  outputs. 
(Pulse  durations  of  20  and  200  times  the  inter-antenna 
propagation  time  are  utilized) . 

Fixed  parameters:   Arrival  Angle  =  60.0° 

SNR  =  Infinite  (no  noise) 
A .   Discussion  and  Conclusions 

The  plots  comprising  Group  1  are  autocorrelation 
and  crosscorrelation  outputs  for  four  decreasing  values  of 
signal-to-noise  ratio.   The  first  series,  shown  in  Figures 
3-2  through  3-11,  is  for  an  ideal  noise-free  signal  environ- 
ment.  The  progression  of  the  matching  coefficient  from 
0.900  to  1.100  in  increments  of  0.05  shows  the  expected 
result  for  the  perfectly  matched  case  (D  =  1.000).   The 
autocorrelation  output  is  an  exact  replica  of  the  input 
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pulse  with  a  gain  factor  of  32,  or  twice  the  array  length. 
The  crosscorrelation  output  is  everywhere  zero  for  the 
matched  case. 

As  the  signal-to-noise  ratio  begins  to  take  on 
finite  values,  however,  the  crosscorrelation  outputs  no 
longer  exhibit  a  zero  value  for  the  matched  case.   It  is 
also  evident  that  the  autocorrelation  outputs  exhibit 
finite  signal-to-noise  ratios  in  the  presence  of  input  noise. 
A  comparison  of  the  matched  autocorrelation  outputs  for 
input  signal-to-noise  ratios  of  10:1,  5:1,  and  2:1  shows 
that  the  output  signal-to-noise  ratios  are  numerically 
equal  to  the  input  signal-to-noise  ratios. 

The  most  important  conclusion  to  be  reached  from 
this  first  group  of  plots,  is  that  the  crosscorrelation 
outputs  still  possess  unique  characteristics  that  distinguish 
them  from  the  outputs  of  the  mismatched  filters,  even  in  the 
presence  of  input  noise.   These  characteristics  are  minimum 
amplitude  and  unipolar  waveforms.   The  peak  amplitudes  of 
the  crosscorrelation  output  waveforms  show  marked  reductions 
as  the  perfectly  matched  case  is  reached.   This  characteris- 
tic is  evident  in  the  plots  in  Group  2  for  array  lengths  of 
10,  8,  and  4  as  well.   In  general,  for  an  input  signal-to- 
noise  ratio  of  5:1,  all  crosscorrelation  outputs  show  an 
approximate  300$  reduction  in  peak  output  as  the  filter 
characteristic  progresses  from  a  mismatched  condition  to  a 
perfectly  matched  condition  (regardless  of  array  length). 
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The  second  identifying  feature  of  the  crosscorrela- 
tion  outputs  is  their  polarity.   For  the  mismatched  case  the 
outputs  exhibit  large  amplitude  noise  peaks  of  both  positive 
and  negative  polarity.   When  the  filters  provide  exact 
matching,  however,  the  output  is  of  a  single  polarity.   For 
the  l6-element  array  utilized  in  Group  1  these  unipolar 
outputs  are  on  the  positive  side.   An  interested  quality  of 
these  unipolar  outputs  was  noticed  during  the  plotting  of 
the  various  outputs  for  this  paper.   That  is,  the  positive 
or  negative  polarity  of  these  perfect  match  crosscorrelation 
outputs  can  be  reversed  by  merely  reversing  the  order  of  the 
complementary  sequences  used  to  phase  code  the  input  signals, 
While  no  attempt  is  made  here  to  justify  this  result,  it  is 
mentioned  because  it  might  be  of  value  in  the  design  of 
recognizers  for  scanning  the  various  correlation  filter 
outputs. 

It  should  also  be  noted  that  the  matched  cross- 
correlation  outputs  do  not  possess  unipolar  qualities  for 
all  array  lengths.   For  instance,  the  matched  crosscorrela- 
tion output  (with  input  noise)  for  an  array  length  of  10, 
as  shown  in  Figure  3-^4,  is  bipolar.   It  was  assumed,  there- 
fore, that  this  characteristic  is  dependent  upon  the  chosen 
code  length.   Not  enough  data  is  available  here  to  attempt 
any  further  explanation  of  this  characteristic. 

The  plots  comprising  Group  2  are  autocorrelation  and 
crosscorrelation  outputs  for  array  lengths  of  10,  8,  and 
4  elements.   These  plots  were  obtained  in  order  to  determine 
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the  effect  of  array  length  on  the  resolution  capability  of 
the  system.   These  plots,  and  those  included  in  Group  1,  can 
be  used  to  derive  an  approximate  maximum  value  of  the 
resolution.   That  is,  for  matching  coefficients  of  0.950 
and  1.050,  which  correspond  to  filters  matched  to  arrival 
angles  1.7°  above  and  below  the  true  arrival  angle,  the 
correlation  outputs  are  indicative  of  a  mismatched  condition. 
Therefore,  the  resolution  of  the  system  is  at  most  1.7°. 
Although  finer  resolutions  are  probable,  and  more  precise 
limits  could  be  assigned  to  this  parameter  by  comparison  of 
additional  output  plots,  it  is  important  to  note  that  the 
resolution  is  not  degraded  by  the  shorter  array  lengths. 
The  basic  conclusion  reached  from  these  results  was  that  a 
system  utilizing  short  (2  or  4  element)  arrays  is  technically 
feasible.   It  would  seem  that  the  only  advantage  to  be  de- 
rived from  using  longer  array  lengths  is  the  increased  signal 
processing  gain  in  the  autocorrelation  mode. 

The  plots  comprising  Group  3  are  autocorrelation  and 
crosscorrelation  outputs  for  pulse  widths  normalized  to 
20  and  200  times  the  maximum  inter-antenna  propagation  time. 
This  propagation  time  is  equivalent  to  the  time  difference 
between  antenna  responses  for  a  wavefront  incident  at  an 
arrival  angle  of  zero  degrees.   It  is  obvious  from  these 
plots  that  the  resolution  of  the  system  is  degraded  for 
signals  of  increased  duration.   The  crosscorrelation  outputs 
which  were  previously  characteristic  of  the  perfectly 
matched  filter,  now  result  for  a  range  of  filter  mismatch. 
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For  20X  pulse  durations  this  output  occurs  over  a  range  of 
filter  mismatches  from  D  =  1.000  (60°  match)  up  to  D  =  1.050, 
which  corresponds  to  an  arrival  angle  of  58.3°.   For  the  200X 
pulse  this  range  of  ambiguous  outputs  runs  from  60°,  for  a 
normally  perfect  match,  up  to  a  matching  coefficient 
corresponding  to  25.8°.   This  last  series  of  output  plots 
indicates  that  the  proposed  system  would  not  be  effective 
against  signals  with  long  duration  modulation  events. 
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AUTOCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  -  16, 
SIGNALAOISE  -  INF.,  PULSE  DURATION  ,-2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITSAlV. 

Figure  3-2 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  -  16, 
SIGNAL/taSE  =  INF.,  PULSE  DURATION  -2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-3 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  -  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT,, 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIVo 

Figure  3-4 
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AUTOCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  -  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  --   16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-5 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-6 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0,900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63,3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/DIV0  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-7 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°)  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSEe  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT,, 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-8 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNALAoISE  =  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  TOUTS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-9 
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CRQSSCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1 .050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (H0RIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/jDIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-10 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/1)IV.  AMPLITUDE  SCALE:   5  UNITS/DIV0 

Figure  3-11 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 


J 


*rT 


:y 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  VJIDTH 100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-12 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-13 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  -  16, 
SIGNAL/NOISE  ■-   10.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-14 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =16, 
SIGNAL/NOISE  -  10.0,  PULSE  DURATION  --   2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.   AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-15 
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AUTOCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME   (HORIZ)    TOUTS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.      AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:      200  UNITS/DIV.     AMPLITUDE  SCALE:      10  UNITS/DIV. 


Figure  3-16 
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CROSSCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY,. 


I 

Jl 

1 

1 

1           n 

!CC 

3C2                        1       p 

U  l| 

J  2  8 

IL 

11 

1. 

TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/t)IV.  AMPLITUDE  SCALE:   $  UNITS/DIV. 

Figure  3-17 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  V/IDTH 100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-18 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-19 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-20 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  ^  16, 
SIGNAL/NOISE  =  10.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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HME  (H0RIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-21 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-22 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°)  ARRAY  LENGTH  =  16, 
SIGNALAOISE  =  5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 


',7  33s!  .33*  330 

TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-23 


101 


AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1 .000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)    TOUTS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.     AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:      200  UNITS/DIV.     AMPLITUDE  SCALE:      10  UNITS/DIV. 


Figure  3-24 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  .-  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV0 


Figure  3-25 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  -  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  -  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-26 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =5.0,  PULSE  DURATION  =  2#0  TLMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-27 
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CROSSCORRELATION  PLOT:  AOA  =  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  16, 
SIGNALAOISE  =  $.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  5  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 

TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-28 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/ftOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


_ 

n 

n 

n 

|}J2 

].■: 

3„G 

::.i 

TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-29 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =5.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (H0RIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-30 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-31 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.900  TIKES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNALAOISE  =  2.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-32 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  -  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  -  2.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-33 
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AUTOCORRELATION  PLOT:      AOA  =  60.0  ,    FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),    ARRAY  LENGTH  =   16, 
SIGNAL/NOISE  =2.0,   PULSE  DURATION   =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)   UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.     AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:      200  UNITS/DIV.     AMPLITUDE  SCALE:      10  UNITS/DIV. 

Figure  3-3-4 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =2.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-35 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  -  16, 
SIGNAL/NOISE  -  2.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-36 
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crosscorrelation  plot:  aoa  =  60.0°,  filter  delay  =  0.900  times  the 
antenna  delay  (filter  is  matched  for  63.3°),  array  length  =  16, 
signal/Noise  =  2.0,  pulse  duration  =  2.0  times  the  antenna  delay. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/blV. 


Figure  3-37 
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CROSSCORRELATION  PLOT:   AOA  --   60.0  ,  FILTER  DELAY  =  0.950  TIKES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  2.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/folV. 

Figure  3-38 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =16, 
SIGNAL/NOISE  ---  2.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (H0RIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-39 
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CROSSCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1 .050.  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =2.0,  PULSE  DURATION  «  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIY.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-40 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  2.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  100  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-4-1 
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CROSSCORRELATTON  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63,3°),  ARRAY  LENGTH  =  10, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 


C   1 

| 

M 
(  ) 
CI 

l ) 

n 

,H 

1 

J-.. 

I 
i 

| 

1  1 

f  > 

C.I 

i 

if- 

i 

ft 

2z: 

i 

1 

1 

Ju   J 

]o6 

330 

C> 

1 

1 

- 

IP 
(  ) 

1 

TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  —  200  UNITS 
EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-42 
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CROSSCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  -  10, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-43 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  -  10, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-44- 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  10, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  TOUTS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-4-5 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  10, 
SIGNAL/frOISE  =5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-4-6 
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CRDSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  8, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE,  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/DIV.  AMPLITUDE  SCALE:  2  UNITS/DIV. 

Figure  3-4-7 
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CROSSCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  8, 
SIGNAL/NOISE  ■-   5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT0 
TIME  SCALE:  200  UNITS/DIV.  AMPLITUDE  SCALE:  2  UNITS/DIV. 

Figure  3-4-8 
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CROSSCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARilAY  LENGTH  -  8, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   2  UNITS/DIV. 

Figure  3-4-9 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  8, 
SIGNAL/NOISE  =5.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSEC  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT 0 

TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   2  UNITS/DIV. 

Figure  3-50 
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CROSSCORRELATION  PLOT:  AOA  =  60.0°,  FILTER  DELAY  =   1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  8, 
SIGNAL/NOISE  -  5.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:  2  UNITS/DIV. 

Figure  3-51 
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CROSSCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  --   0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  4, 
SIGNAL/NOISE  =5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIVo  AMPLITUDE  SCALE:   1  UNIT/DIY. 

Figure  3-52 
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CROSSCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  /,, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  -  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   1  UNIT/DIV. 

Figure  3-53 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  U, 
SIGNAL/NOISE  -  5.0,  PULSE  DURATION  =  2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/TjIV.  AMPLITUDE  SCALE:   1  UNIT/DIV. 

Figure  3-54- 
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CROSSCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  -  4, 
SIGNAL/NOISE  =  5.0,  PULSE  DURATION  --   2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  200  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   1  UNIT/DIV0 

Figure  3-55 
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CROSSCORRELATION  PLOT:   AOA  -  60.0°,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  4, 
SIGNAL/NOISE  =5.0,  PULSE  DURATION  =2.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  - —  200  UNITS 
EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   1  UNIT/DIV. 

Figure  3-56 
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AUTOCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  =  0.600  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  72.6°),  ARRAY  LENGTH  =16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 


Juu 


]C2 


Ju'i 


use 


]GlJ 


TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 4.00  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-57 
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AUTOCORRELATION  PLOT:  AOA  =  60.0°,  FILTER  DELAY  =  0.800  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  66.4°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-58 


136 


AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAI  -  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV« 


Figure  3-59 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (H0RIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 


Figure  3-60 
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AUTOCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  -  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT „ 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-61 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  -  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  ^00  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:  200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-62 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1.100  TIMES  THE 
ANTENNA  DELAI  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  ---  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  ^00  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-63 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.200  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  53.1°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  —  400  UNITS 
EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-64- 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1 .400  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  45.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-65 
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CR0SSC0RRELAT10N  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.900  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  63.3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-66 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.950  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  61.7°),  ARRAY  LENGTH  -  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-67 
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CROSSCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  _  400  UNITS 
EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 

Figure  3-68 
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CROSSCORRELATIOH  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.050  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  58.3°),  ARRAY  LENGTH  -  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIVo 

Figure  3-69 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =1.100  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  56.6°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  20.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   5  UNITS/DIV. 


Figure  3-70 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.200  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  84,3°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  VJIDTH  400  UNITS 

EQUAL  ONE  PULSEC   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-71 
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AUTOCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.600  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  72.6°),  ARRAY  LENGTH  r.  16, 
SIGNAL/NOISE  -  INF.,  PULSE  DURATION  --  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 400  UMTS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-72 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-73 
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AUTOCORRELATION  PLOT:   AOA  -  60.0  ,  FILTER  DELAY  =  1.400  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  45.6°),  ARRAY  LENGTH  -  16, 
SIGNAL/NOISE  --   INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  £00  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DTV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-74 
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AUTOCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  --   1.800  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  25.8°),  ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-75 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  0.200  TIMES  TIE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  &U3°),  t^^1   LENGTH  -  16, 
SIGNAL/NOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV0  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-76 
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CROSSCORRELATION  PLOT:   AOA  =  60.0  ,  FILTER  DELAY  =  0.600  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  72.6°),  ARRAY  LENGTH  =16, 
SIGNALAOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHTo 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-77 
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CROSSCORRELATION  PLOT:      AOA  -  60.0  ,    FILTER  DELAY  =  1.000  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  60.0°),    ARRAY  LENGTH  =  16, 
SIGNAL/NOISE  =  INF.,    PULSE  DURATION  =  200.0  TIRES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)   UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 4-00  UNITS 

EQUAL  ONE  PULSE.     AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:      200  UNITS/DIV.     AMPLITUDE  SCALE:      10  UNITS/DIV. 

Figure  3-78 
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CROSSCORRELATION  PLOT:      AOA  =  60.0°,    FILTER  DELAY  =  1.400  TIMES  THE 
ANTENNA  DELAY   (FILTER  IS  MATCHED  FOR  4.5.6°),    ARRAY  LENGTH  =   16, 
SIGNAL/NOISE   =  INF.,    PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH  ^00  UNITS 

EQUAL  ONE  PULSE.   AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-79 
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CROSSCORRELATION  PLOT:   AOA  =  60.0°,  FILTER  DELAY  =  1.800  TIMES  THE 
ANTENNA  DELAY  (FILTER  IS  MATCHED  FOR  25.8°) ,  ARRAY  LENGTH  =  16, 
SIGNALATOISE  =  INF.,  PULSE  DURATION  =  200.0  TIMES  THE  ANTENNA  DELAY. 
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TIME  (HORIZ)  UNITS  ARE  NORMALIZED  TO  INPUT  PULSE  WIDTH 400  UNITS 

EQUAL  ONE  PULSE.  AMPLITUDE  UNITS  ARE  NORMALIZED  TO  PULSE  HEIGHT. 
TIME  SCALE:   200  UNITS/DIV.  AMPLITUDE  SCALE:   10  UNITS/DIV. 

Figure  3-80 
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IV.   SUMMARY  AND  RECOMMENDATIONS 

A.   SUMMARY 

The  phase  coding  radio  direction  finder  using  comple- 
mentary sequence  pairs  has  been  evaluated  with  the  objective 
of  applying  the  basic  technique  to  a  realizable  model  of  the 
system.   A  sinusoidal  analysis  of  the  carrier  phase  relation- 
ships which  exist  in  the  various  portions  of  the  signal 
process  has  shown  that  the  unique  autocorrelation  and  cross- 
correlation  outputs  will  be  produced  for  any  arrival  angle, 
if  the  system  is  matched  for  that  angle. 

It  was  shown  that  the  inherent  ambiguities  of  a  single 
array  system,  caused  by  elevation  angle  and  adjacent  quadrant 
azimuth  components,  can  be  resolved  by  using  a  common  crossed- 
array  interferometer  technique.   This  technique  requires  that 
a  set  of  correlation  filters  corresponding  to  a  180°  range 
of  arrival  angles  be  available  for  each  array. 

A  computer  model  of  the  system  was  used  to  evaluate  the 
performance  of  the  signal  processor  under  varying  input  noise 
conditions,  and  for  varying  array  lengths  and  signal  duration 
parameters.   The  conclusions  reached  from  this  analysis 
indicate  that  the  angular  resolution  of  the  system  will  not 
be  seriously  affected  by  the  array  length  or  input  signal-to- 
noise  ratio.   That  is,  for  most  input  signals,  the  resolving 
capability  of  the  system  is  dependent  only  upon  the  number 
and  incremental  differences  between  the  time  delay  segments 
of  the  correlation  filters  used.   However,  it  was  also 
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determined  that  the  duration  of  the  signal  model  used  for 
this  analysis  had  a  very  significant  effect  upon  the  resolu- 
tion. For  signal  pulse  widths  that  were  large  compared  with 
the  time  of  arrival  difference  between  antennas,  the  resolu- 
tion of  the  system  was  seriously  degraded.  As  a  result,  it 
was  concluded  that  the  system  would  not  be  effective  against 
signals  with  this  characteristic. 

B.   RECOMMENDATIONS  FOR  FURTHER  INVESTIGATION 

Since  the  analysis  of  the  proposed  system  was  based 
entirely  upon  the  relative  time  of  arrival  of  distinct 
modulation  events,  and  since  this  was  shown  to  place  a 
serious  limitation  on  the  type  of  modulation  signals  that 
can  be  processed  by  the  system,  an  alternate  technique 
should  be  pursued.   That  is,  the  utilization  of  carrier 
phase  differences,  resulting  from  the  propagation  delay 
between  antennas,  might  be  a  more  acceptable  parameter  to 
use  rather  than  time  difference. 

A  second  recommendation  would  be  the  construction  of  a 
simple  scale  model,  utilizing  one  or  two  small  arrays  in 
order  to  evaluate  the  operation  of  the  system  in  a  multi- 
signal  environment. 
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APPENDIX 
COMPUTER  PROGRAM 

The  computer  simulation  is  a  digital  program  written  in 
FORTRAN  IV  level  G.   It  was  processed  by  an  IBM  360/67 
system  and  the  built  in  functions  are  those  normally 
available  on  this  system. 

The  program  requires  a  data  deck  input  for  the  specific 
complementary  code  used  in  addition  to  the  parameters 
required  to  process  the  code.   The  following  words  are 
peculiar  to  this  program. 

C0DE1:   First  sequence  In  a  complementary  code 

C0DE2 :   Second  sequence  in  a  complementary  code 

SIGN:   Summation  sign,  used  for  addition  of  correlator 
outputs 

DELDEG:   Incremental  angle  shift  in  degrees,  used  to 
obtain  a  new  value  for  the  angle 

D:   Filter  delay  starting  value 

START:   The  starting  value  in  degrees  of  the  incoming 
signal 

STOP:   The  final  value  in  degrees  of  the  incoming  signal 

START-T:   The  starting  time  for  the  computation 

STOP-T:   The  completion  time  for  the  computation 

DSTEP:   Incremental  step  for  the  parameter  D,  used  to 
obtain  a  new  value  for  the  filter  (correlator) 
delay 

NRANTS:   The  number  of  antennas  in  use 

SPACE:   Spacing  between  antennas;  antennas  have  uniform 
spacing 

DELT:   Computation  time  increment 
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TIMES:   Multiplying  factor  used  in  the  subprograms  in 
order  to  utilize  integer  arithmetic  and  logic 

J J :   Control  for  suppressing  data  printout;  do  not 
print  if  JJ^2 

JJJ:   Control  for  correlation  function;  if  JJJ=1  then 
cross-correlate ,  if  not  then  autocorrelate 

JJJJ:  Control  to  suppress  the  final  plot  which  is  a 
summary  of  maximum  amplitude  versus  angle;  if 
JJJJ=2  then  do  not  plot 

The  data  deck  must  be  assembled  in  a  specific  sequence 
in  order  to  have  the  program  operate  properly.   Additionally 
if  the  data  is  not  known  a  card  must  still  be  loaded  into 
the  data  deck  in  order  not  to  invalidate  the  read  instruc- 
tions.  The  data  must  be  provided  whether  valid  or  not. 
The  data  cards  are  assembled  as  in  Figure  A-l.   The  first 
card  in  the  data  deck  is  C0DE1  and  the  last  is  JJJJ.   All 
data  cards  are  punched  starting  in  column  one.   C0DE1  and 
C0DE2  are  formatted  floating  point  F4 . 1  and  each  card  can 
have  up  to  eighteen  entries.   Each  code  is  read  from  one 
card,  if  longer  codes  (greater  than  18)  are  required  then 
the  read  statement  must  be  changed  to  read  more  than  one 
card  for  each  code.   NRANTS,  JJ,  JJJ  and  JJJJ  are  formatted 
in  integer  111  and  all  other  data  cards  are  formatted  in 
floating  point  P11.6. 

The  program  can  provide  for  a  pulse  signal  input  with 
the  pulse  duration  equal  to  any  multiple  of  the  delay  due 
to  antenna  separation.   The  program  will  automatically  scale 
the  built  in  pulse  duration  (W)  and  the  antenna  separation 
(SPACE)  to  accommodate  any  code  length.   Additionally  the 
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computation  time  is  automatically  scaled  in  order  to 
include  all  the  correlation  outputs  greater  than  zero  and 
also  provide  for  the  maximum  number  of  computations  per 
angle  and  filter  delay  setting.   The  maximum  number  of 
computation  intervals  for  a  fixed  angle  and  delay  is  2000. 
The  codes  are  read  into  computer  memory  and  all  necessary 
operations  are  accomplished  in  order  to  provide  the  proper 
autocorrelation  and  crosscorrelation  functions  as  a  function 
of  time,  angle  and  filter  delay.   If  data  Is  provided  which 
is  beyond  the  array  storage  set  aside  for  the  computation 
then  the  program  will  attempt  adjustments  in  the  size  of  the 
data.   After  twenty  attempts  with  no  success  the  program  will 
terminate  and  print  out:   DATA  EXCEEDS  ARRAY  STORAGE... 
JOB  TERMINATED. 

The  program  will  input  all  parameters,  compute  other 
required  parameters  in  accordance  with  the  inputted  param- 
eters, compute  an  output  for  a  fixed  time,  angle  and  delay 
and  then  step  the  time  variable  to  the  next  increment.   After 
the  time  variable  has  been  stepped  across  its  range  the 
filter  matching  coefficient  (D)  is  incremented  to  a  new 
value,  and  the  output  is  computed  for  all  time  values  as  the 
time  variable  is  again  swept  through  its  range.   After  the 
matching  coefficient  has  been  swept  through  values  from  its 
initial  value  to  a  final  value  equally  spaced  above  D=1.000, 
the  angle  variable  may  then  be  incremented  to  a  new  value. 
The  outputs  are  then  computed  across  the  time  range  and  the 
angle  range.   In  this  manner  the  outputs  of  the  filters 
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(correlator  outputs)  are  computed  for  all  times  of  interest, 
all  angles  of  interest  and  all  desired  variations  of  the 
filter  delay  lines. 

The  outputs  available  from  this  program  are  printer 
plots  and  fine  grain  continuous  plots  of  correlation  output 
versus  time  for  each  and  every  value  of  angle  and  filter 
delay.   The  printer  plots  contain  all  pertinent  parameters 
and  scaling  information. 

A  listing  of  important  parameters  and  statements  con- 
cerning the  interpretation  of  the  various  control  commands 
are  printed  at  the  beginning  of  each  computer  run  to  aid  in 
analysis  of  the  data  and  the  plots.   The  maximum  core  size 
is  150K  and,  due  to  the  fine  grain  computation  required  by 
the  program,  execution  time  for  most  investigations  is  in 
the  15  to  20  minute  region. 
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TYPICAL  DATA  DECK 


1.0  1.0  1.0-1 

1.0  1.0  1.0-1 

1.000000 

3.000000 

0.990000 

60.000000 

60.000000 

0.000000 

0.000000 

16 

2 

2 

2 


1.0 
1.0 


0-1 
0-1 


1.0  1.0  1 
1.0-1.0-1 


1.0-1.0-1 
■1.0  1.0  1 


■1.0  1 

1.0-1 


0-1.0 
0  1.0 


CARD  1. . .FIRST  SEQUENCE  OF  THE  CODE  PAIR (CODE1) 

CARD  2... SECOND  SEQUENCE  OF  THE  CODE  PAIR (CODE2) 

CARD  3-  •  .SUMMATION  SIGN (SIGN) 

CARD  l\.  .  .ANGLE  STEP (DELDEG) 

CARD  5.  .  .FILTER  DELAY  STARTING  VALUE (D) 

CARD  '6.  .  .ANGLE  START (START) 

CARD  7...  ANGLE  STOP (STOP) 

CARD  8.  .  .COMPUTATION  TIME  START (START-T) 

CARD  9.  .  .FILTER  DELAY  STEP (DSTEP) 

CARD  10  .  .  NUMBER  OF  ANTENNAS (NRANTS  ) 

CARD  11..  CONTROL  FOR  DATA  PRINT (JJ) 

CARD  12.. CONTROL  FOR  TYPE  OF  CORRELATION (JJJ) 

CARD  13.. CONTROL  FOR  OUTPUT  VS  ANGLE  PLOT (JJJJ) 


FIGURE  A-l 
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